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ABSTRACT 


Recently, a new method based on Feed Forward Neural Networks (FFNNs), 
called the 'Delta' method, has been proposed for estimating the aircraft stability and 
control derivatives (parameters) from flight data In the past, the Delta method has been 
validated on simulated data and a few limited set of real flight data only for lateral- 
directional motion. In the present work, the Delta method has been validated on 
extensive real flight test data, both for longitudinal and lateral-directional motion. The 
Delta method has been firstly used to estimate parameters from raw flight data that may 
contain bias, scale factor and measurement noise errors. Further, the Delta method has 
also been used to extract parameters from the corrected flight data (that is flight data 
obtained after applying data con^atibility checks to raw flight data). The results of the 
Delta method show good consistency and low standard deviation. Further, the estimated 
aerodynamic coefficients and estimated response obtained from the estimates, compare 
very well with the respective measured quantities. Thus, the Delta method estimates are 
feirly accurate. The results show that the Delta method can be advantageously applied to 
estimate parameters from real flight data. The Delta method, unlike conventional 
methods, neither requires a priori fix of model nor involves integration of equations of 
motion of the system nor requires initial guess values of the aerodynamic parameters to 
be given. This is a distinct advantage of the Delta method. Further, the results show that 
the Delta method is enable of estimating parameters from the raw flight data using the 
measured values of the motion variables that are not corrected for bias and scale factor 
errors. 
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CHAPTER 1 


INTRODUCTION 


Airplane parameter estimation is the process of estimating numerical values of the 
aerodynamic stability and control derivatives from flight test data. Accurate values of 
these parameters are required in the mathematical model for description of the airplane 
dynamics. Such models are useful for specification of stability augmentation systems, 
flight simulators and fi)r verification of theoretically predicted models. 

There are three distinct approaches for estimating stability and control derivatives: 

• Theoretical Methods 

• Wind tunnel testing 

• Flight Testing 

At the preliminary design stage of an airplane, theoretical method^*^’^ are useful in spite 
of their limited accuracy. The wind tunnel test improves the accuracy of estiniation but it 
is a time consuming and expensive way of estimating stability and control derivatives. 
Precise simulation of control surfaces, power effects and flight conditions is difficult. 
The model tested in wind tunnel is generally slightly different from the actual airplane 
due to last minute configuration changes. Other reasons for discrepancies between flight 
and wind tunnel results are Reynolds number discrepancies and interference due to 
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support system. It is therefore desirable that the wind tunnel estimates be corroborated 
with the estimates from the actual flight test data. 

The most commonly used methods for parameter estimation from flight data are 
broadly classified into the following categories. 

• Equation error method 

• Output error method (the maximum likelihood method being the most popular) 

The principle of least squares is used in the equation error method wherein the error 
gets minimized with respect to the unknown parameter in each of the equations'^. Its 
advantages include computation simplicity, non-iterative nature and applicability to both 
linear and non-linear models. In the output error method, the error between the measured 
and the model response produced by identical inputs is minimized. The method assumes 
that there exist no modeling errors. This method processes the measurement noise while 
assuming the model representation of the given system to be exact. The maximxim 
likelihood (ML) estimator and its several variants have been the most successfully used 
methods for estimating parameters from flight data^’^. However, ML method and its 
variants need a priori fix of an aerodynamic model. Further, for highly maneuverable 
aircraft (unstable in open loop), ML method feces difficulties during integration of 
system equations and leads to diverging solutions^. 
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The above mentioned methods require accurate measurement of input output 
variables. This requirement is difficult to meet in practice due to many types of 
measurement errors in flight data*’^. 

A few typical measurement errors present in data are: the bias (also called zero 
shift), scale factor, time shift and measurement noise. For conventional parameter 
estimation methods, it is, therefore, required that the data compatibility check (sometimes 
referred to as flight path reconstruction) be performed to get accurate results. Thus, the 
data compatibility check becomes an integral part of the process of parameter estimation. 

The data compatibility is performed using the standard kinematic equations, and 
applying a ML algorithm®. The data compatibility check provides accurate information 
about the aircraft states, and also estimation of biases, scale fectors and time shifts errors 
in the recorded data. Elimination of such errors from the flight measurements prior to 
estimation of parameters helps to improve the accuracy of the estimates. However, the 
process of data compatibility check for estimation of errors can itself be quite sensitive to 
the methodology employed, and at times, it may be as time consuming as the process of 
estimating parameters itself’*’®. 

More recently, investigations have been carried out to ejq)lore the potential of 
Artificial Neural Networks (ANNs) for aircraft aerodynamic modeling and parameter 
estimation. ANNs are biologically inspired, massively parallel conq)utational paradigms 
that mimic, to some degree, actual neurological systems. Historically, ANNs were 
developed out of continued interest in the constructional and fimctional details of the 
human brain, and the research there upon to understand its neurological capabilities such 
as its ability to learn, to retain and to recall information on demand. The ANNs have 
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been successfully used for modeling in various fields such as signal processing, pattern 
recognition, system identification and control, and so on. 

A class of artificial neural networks called the ‘feed forward neural networks’ 
(FFNNs) are most widely used among all types of ANNs for aircraft modeling. As plenty 
of literature is available on FFNNs, here we present only a brief outline about them. The 
FFNNs have neurons arranged in layers like directed graphs, and these are static in 
nature. FFNNs lead to a black box model in which no physical significance can be 
assigned to either the network structure or to the network weights. FFNNs work as 
general function approximators, and are capable of approximating any continuous 
function to any desired accuracy by an appropriate network architecture^®. A brief write 
up on the various aspects of FFNNs is presented in Appendix 'A'. 

The ability of FFNNs has been utilized to model aircraft aerodynamics 
successfully^ Very recently, Raisinghani et al'^’*® have used FFNNs for modeling 
aircraft aerodynamics wherein the measured responses and control inputs are mapped to 
the aerodynamic coefficients; such a neural model is subsequently used to estimate 
aircraft parameters by the Delta method*^’^^’** *’. It is also mentioned that the FFNN 
modeling does not require any a priori postulation of the mathematical model for the 
aircraft aerodynamics. Applicability of the Delta method for estimation of lateral- 
directional parameters from simulated as well as limited real flight data has been 
demonstrated, respectively, by Ghosh et al*’’^* and Raisinghani et al^^’^®. Further, in 
Ref. 20, it was shown that the Delta method can be applied advantageously for extracting 
aircraft parameters from flight data containing unknown bias, scale factor and 
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measurement noise errors. The distinct advantage of using the Delta method as reported 
in the above reference, can be summarized as follows. 

(a) . The Delta method can be regarded as a non-parametric modeling method; both the 

structure and parameters need not be known a priori. 

(b) . The Delta method can be applied directly for extracting aerodynamic parameters 
from flight data containing unknown bias, scale factor and measurement noise errors. 

(c) For the Delta method, there is no need to estimate the initial conditions of flight 
variables. 

Raisinghani et al^^’*^ have demonstrated the above advantages of the Delta 
method by applying it to a limited sets of real flight data (only for lateral-directional 
motion) which were obtained from abroad. The next obvious step was to further 
strengthen the efficacy and accuracy of these methods by validation on actual flight data 
for longitudinal motion. The Delta method was also required to be re-validated on 
sufficient sets of real flight data generated for lateral-directional motion by giving 
specific test inputs. 

Thus file motivation for the present work was to carry out a well planned, 
conqirehensive and exhaustive flight trial exercise to gather sufficient real flight test data 
for validating the Delta method with respect to the advantages mentioned in paragraphs 
(a), (b) and (c) above. Towards this end, a joint project was undertaken by IIT Kanpur 
and Aircraft and Systems on Testing Establishment (ASTE), Indian Airforce (TAF), 
Bangalore, such that the flight trials would be undertaken at ASTE and subsequent 
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analysis would take place at IIT Kanpur. Accordingly, Kiran MK lA aircraft, available 
with ASTE was chosen as the test aircraft. 

It must be noted that such an exercise was undertaken for the first time in the 
country, wherein, we have designed our own experiments to gather exhaustive real flight 
test data for parameter estimation using the Delta method. 

The accuracy and reliability of parameter estimates obtained via any of the 
estimation methods, depends heavily on the amoimt of information available in the flight 
data. In the seventies, Mehra, Gupta, et have significantly contributed to the 

optimal input design. Although sporadic developments on optimal input design have 
been reported since then, the 3-2- 1-1 input (explained in subsequent chapter) remains to 
date as the one most accepted worldwide and utilized by the international flight test 
co mmunit y. 

Accordingly, the majority of the flight test data on the Kiran aircraft was 
generated using the 3-2-1-1 control input. Apart fi'om this, other forms of inputs, such as 
doublet and pulse inputs, were also given to enable a study of effect of input form on 
parameter estimates. The flight data thus generated will necessarily have measurement 
errors like bias and scale fectors etc. These sets of data are referred to as raw flight data 
for future reference. The next step was to correct the flight data by incorporating the 
appropriate values of bias and scale factors obtained though data coir 5 )atibiiity check and 
these will be referred to as 'corrected data' in the thesis. 

In the present work, the Delta method is used for extracting parameters from real 
flight data corresponding to longitudinal and lateral-directional dynamics. The Delta 
method is used for extracting parameters Ifrom raw as well as corrected flight data. The 
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results show that the Delta method can be advantageously applied to estimate parameters 
from real flight data. This method neither requires a priori fix of model nor involves 
integration of equations of motion of the system. Further, the results show that the Delta 
method is capable of estimating parameters from the raw flight data using the measured 
values of the motion variables that are not corrected for the bias and/or the scale factor 
errors. 

The work carried out in this thesis is presented in four chapters. In Chapter 1, a 
brief introduction is given on the subject of aircraft parameter estimation. Additionally, 
the motivation for the present work is also discussed. Chapter 2 deals with the 
maneuvers used for generation of simulated and real flight data. Further, the details of 
the data con:q)atibility checks and the Delta method of FFNN are outlined. Chapter 3 
consists of the descriptive study of the results obtained. Finally, Chapter 4 gives the 
conclusions drawn from the present work and also the recommendations for future work. 



CHAPTER 2 


GENERATION OF FLIGHT DATA AND DATA COMPATIBILITY CHECK 

2.1 GENERAL 

This chapter describes the various types of simulated and real flight data 
generated for use in aerodynamic modeling and parameter estimation via the Delta 
method using FFNN. The procedure used for generating the simulated data is &st 
ejqplained in brief Next, the work undertaken to generate real flight data in coordination 
with Indian Air Force (lAF) is described. The preprocessing procedure used for the real 
flight data is also discussed. 

The symmetry about the X-Z plane of the aircraft and small angle approximations 
allow separation of the equations of motion into two largely independent sets describing 
the longitudinal and the lateral-directional motion of aircraft. The imknown parameters 
also separate into longitudinal and lateral-directional sets. We have assumed the 
separation of longitudinal and the lateral-directional dynamics in the present work. 
Accordingly, simulated data generation, aerodynamic modeling and parameter estimation 
are carried out separately for the longitudinal and the lateral-directional motions. For 
generating the simulated flight data, the equations of motion pertaining to longitudinal 
and lateral-directional dynamics were solved using fomth order Runga-Kutta method. At 
this stage, it must be emphasized that we need neither the equations of motion nor their 
solutions for the FFNN modeling and subsequent parameter estimation via the Delta 
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method . The procedure used for generation of simulated data is now discussed in the 
next paragraph. 


2.2 SIMULATED FLIGHT DATA -LONGITUDINAL MOTION 

As a first step, to check the proper working of all the FFNN programs written in 
MATLAB, it was decided to use simulated flight data initially. This would also increase 
the confidence level before working on real flight data. It would also provide a basis for 
comparison of the results as well as the effect of various variables when dealing with real 
flight data. The simulated flight data was obtained by solving following equations of 
motion corresponding to short period mode (u = constant)^'*. 


d - q = -[puS/2m]CL 

(2.1) 

q = [pu^Sc/2Iy]Cni 

(2.2) 


where. 


Cl — Clo + CLa<^ + CLq + C[^de 


Cm ~ Cjho + + Cjjiq + CjjjSgSe 

For generation of simulated flight data, the Kiran 
example aircraft. The flight conditions, mass, 
characteristics as given in Table 2.1 below were fed 


MK LA aircraft was chosen as the 
moment of inertia and geometric 
into the above equations along with 


the aerodynamic parameters. 


if 


9 



Table 2.1 Flight Conditions, Mass, Moment of Inertia and Geometric Characteristics of 
Kiran Aircraft 


True Air Speed 

V = 

103.2157 m/sec 

Mass 

m = 

3442.5 Kg 

Density 

p= 

0.7 Kg/m^ 

Wing Area 

s = 

19 m^ 

Mean Aerodynamic chord 

C = 

1.891 m 

Wing span 

b = 

10.7 m 

Half wing span 

s=b/2 = 

5.35 m 

Moment of Inertia 

Ix = 

4332.084 Kg-m^ 


Iy = 

15132.802 Kg-m^ 


Iz = 

18548.354 Kg-m^ 


Ixy = 

20.517 Kg-m^ 


The longitudinal aerodynamic parameters were taken from Ref. 25, and are as 
shown in Table 2.2 below. These values will be henceforth referred to as true values for 
the longitudinal stability and control derivatives. The values of the parameters Cto and 
Cmo were not available and hence they were assumed to be zero. 


Table 2.2 Aerodynamic characteristic of BCiran MK 1 A Aircraft - Longitudinal Motion 


SL 

No. 

Parameter 

Value used (per radian) 

1. 

CLa 

4.9 


Clq 



ClSc 

0.4 


Cma 

-0.60 

5. 

Cmq 

-13.0 

6. 

Cmde 

-0.68 


The control ii^ut form used to generate the simulated data for longitudinal motion 
was the flight 3-2-1-1 Se input with an amplitude of 0.1 radian and time period of seven 
seconds as given in Fig. 2.1 a 
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Rudder Defln. (radian) 



Fig. 2.1a Plot of Flight 3-2-1 -1 Elevator Input Used For 
Generation Of Simulated Data 



Fig. 2. 1 b Plot of Flight Rudder Doublet Input Used For 
Generation Of Simulated Data 
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2.3 SIMULATED FLIGHT DATA - LATERAL DIRECTIONAL MOTION 


The simplified equations of motion used for generating the lateral-directional 


simulated flight data are as follows^*: 


P = — r + (g/V)sin(j) + (qS/mV)Cy 

(2.3) 

p = {l/(I^I^ - I^xz)}{qSs(I,C, + I^C„)) 

(2.4) 

r = {1/aA -I^xz)}{^s(I^C„ +I^C,)} 

(2.5) 

<i>=p 

(2.6) 

0 = — rsin <f> 

(2.7) 


where (3 = vPV 

q = (l/2)p(u^ -1- + w^) 

V= (u^ + v^+w^y^^ 

c I = c,„ + C,p(^) + Ct (^) + Cipp + CeaSa + C b A 

— Cjjo + Cnp(— ) + Cjj]-(— ) + Cjipp + Cjjgada + Cjjgfdr 

Cy = Cyo + Cyp(— ) + Cy|-(— ) + Cy^P + Cy5a8a + Cyg|-5r 

In this case, the Kiran MK lA was chosen as the example aircraft and hence the same 
mass/geometric characteristic as given in Table 2.1 were used in the Eq. (2.3) to (2.7). It 
may be mentioned that lateral-directional aerodynamic parameters of BCiran aircraft were 
not readily available. Therefore, to generate simulated flight data, tentative values based 
on analogous comparison with similar sized aircraft were used. These tentative values as 
shown in Table 2.3 below will be henceforth called as the true values for the lateral- 
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directional stability and control derivatives. In this case, since the parameters Cio , Cno 
and Cyo were unknown, they were assumed to be zero. 


Table 2.3 Aerodynamic characteristic of the Kiran aircraft - Lateral-Directional Motion 


SL No. 

Parameter 

Value used (per radian) 

1. 

Cto 

-1.203 

2. 

C,r 

0.189 

3. 

Ci3 

-0.144 

4. 

C,8r 

0.032 

5. 

Cno 

-0.210 

6. 

C„r 

-0.805 

7. 

Cnj3 

0.193 

8. 

Cn5r 

-0.165 

9. 

Cyp 

0.990 

10. 

Cyr 

1.594 

11. 

CvB 

-0.773 

12. 

Cydr 

0.302 


The control ii^ut form used for generating the simulated flight data for lateral-directional 
motion was the flight 6r doublet input as presented in Fig. 2.1b. Having described the 
procedure for the generation of simulated flight data, the next paragraph gives the details 
of the procedure followed to generate real flight data by performing different types of 
maneuvers using Kiran aircraft. 

2.4 GENERATION OF REAL FLIGHT DATA 

Aircraft and Systems Testing Establishment (ASTE) of lAF is a premier flight test 
agency in the country. It is responsible for undertaking all flight test activities of Indian 
Armed Forces. This unit located at Bangalore, is also additionally responsible for 
training of test pilots and flight test engineers. Hence, it was decided to conduct the flight 
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trials for the present work at ASTE, lAF. The test aircraft chosen for the flight trials was 
the lAF trainer aircraft called as 'KIRAN MARK lA'. This aircraft is a twin seater low 
wing aircraft fitted with a single jet engine. 

The standard instrumentation scheme presently available with ASTE, lAF was 
used for flight data acquisition for this work. The list of parameters acquired are as listed 
in Table 2.4 below. The instrumentation scheme is a digital data acquisition system 
based on pulse code modulation (PCM) concept. The resolution and accuracy of the 
instrumentation scheme is as shown in Table 2.4. 


Table 2.4 List of Parameters Acquired by the PCM based DAQ Instrumentation System 



Parameter 

Symbol 

Unit 

Resolution 

Sampling 

Rate 

1. 

Altitude 

Zp 

millibar 

0.256 

16 

2. 

IAS 

Vi 

miUibar 

0.097 

16 

3. 

AOA 

a 

degree 

0.034 

16 

4. 

AOS 


degree 


16 

5. 

Elevator Position * 

5e 

degree 


16 

6. 

Aileron Position * 

6a 

degree 

0.012 

16 

7. 

Rudder Position * 

5r 

degree 


16 

8. 

Pitch Rate 

q 



16 

9. 

Yaw Rate 

r 



16 

10. 

Roll Rate 

p 



16 

11. 

Longitudinal Load Factor 

nx 

g 


16 

12. 

Lateral Load Factor 

% 

g 


16 

13. 

Normal Load Factor 

Dz 

g 

0.0023 

16 

* Control 

Sur&ce Positions 


A total of five sorties were flown each of approximately one hour duration in 
calm and fine weather conditions. The details of the sorties are as shown in Table 2.5 . 
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Table 2.5 Details of the Sorties Flown for Generating Real Flight data 



Date Flown 

Duration 

Purpose 

FLT 1 

06“ July 2000 

01:05h 

Shake down sortie for the instrumentation 
scheme. 

FLT2 

06*^ July 2000 

OLOOh 

Angle of attack calibration in flight 

FLT 3 

07“ July 2000 

OLlOh 

To excite the longitudinal response of the 
aircraft 

FLT 4 

10^^ July 2000 

00:55h 

To excite longitudinal response of the 
aircraft. 

FLT 5 

n*** July 2000 

OLlOh 

To excite lateral-directional response of 
the aircraft. 


The importance of choosing appropriate inputs (control surface deflections) for 
exciting specific modes of an aircraft or executing specific maneuvers during flight 
testing for parameter estimation has been well recognized^. The accuracy and reliability 
of parameter estimates depends heavily on the amount of information available in the 
aircraft response. It is necessary that the correct mode of aircraft response is excited by 
giving carefully selected inputs. 

Based on Ref 9, the following test inputs were used to generate the aircraft flight 

data:- 


(a) Pulse input of different amplitudes and time period. 

(b) Doublet inputs. 

(c) 3-2- 1-1 inputs. 

(d) Roll maneuvers as described in Ref 9. 

Although, it is well known that the 3-2-1-1 input is the most optimal input,^*’^’^^ 
the other input forms were used so as to initiate a study on sensitivity of parameter 
estimates to different input forms. However, due to paucity of tune, this aspect could 
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not be studied in detail. Hence in the present work the flight data generated using the 
3-2-1 -1 elevator input for longitudinal motion and doublet rudder/aileron input for lateral 
directional motion have been used for FFNN modeling and subsequent parameter 
estimation by the Delta method. The data sets which could not be analyzed in this thesis 
will be taken up for study in future. 

The test altitude was chosen as Zp =12,000 feet with a tolerance band of ±200 feet 
for all the maneuvers. The speeds were chosen so as to cover the cruise range of the 
aircraft, namely 160, 180, 200 and 220 knots, indicated airspeed (KIAS). All the 
maneuvers were undertaken for undercarriage retracted and no flap deflection 
configuration. 

For the present work, the data sets obtained, pertaining to different maneuvers or 
the test points executed, are specified by the indicated airspeed (Vj). Because more than 
one maneuver or test point were done/performed at the same speed, they are designated 
additionally by the flight number and the sequence number in which the maneuvers were 
executed. Accordingly, a simple nomenclature system is formulated to refer to the 
various data sets analyzed. Each set of data is referred to by the flight number, 
designated as FLT x where x=3, 4 or 5 followed by a slash and then the corresponding 
test point number (that is, sequence number of the maneuver) written as TPy where y = 1, 
2, 3 etc. Thus FLT 4 /TP 2 refers to the data set obtained in flight number 4, test point (or 
maneuver) number 2. 

Tlie real flight data thus generated was then analyzed. The first step in the 
analysis was to carry out the data con:q)atibility check "wflich is described in the next 
paragraph. 
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2.5 DATA COMPATIBILITY CHECK 


Real flight data generally have unknown biases, scale fectors and time shifts in 
the recorded data^’'*. Bias or zero shift errors take place due to drift in the transducer 
from its zero position after calibration. This will be a constant value error. The scale 
error takes place due to an interference effect during the measurement process. This 
causes a constant scale error to be introduced in the output of the transducer. These 
errors cause data incompatibility. Data incompatibility could arise by way of the 
measured incidence angles not being in agreement with those reconstructed from the 
accelerometer and rate gyro measurements. Therefore, it is essential that the data 
compatibility check be carried out before using the data for aerodynamic modeling and 
parameter estimation. In other words, the data compatibility check, which is also called 
as the flight path reconstruction, is an integral part of aircraft parameter estimation. For 
conventional methods of aircraft parameter estimation, it is well known that data 
compatibility checks prior to estimation of derivatives helps to improve the accuracy of 
the estimates’. Using the well defined kinematic model, the data con^)atibility check 
provides an accurate estimates of the aircraft states, and estimates of biases, scale fectors 
and time shifts in recorded data ’ . 

To estimate the aircraft states, the unknown scale &ctors and biases of the state 
variables for the real flight data of Kiran aircraft, a suitable 'Maximum Likelihood 
Estimation' program was developed (in Matlab 5.3). In applying the MLE technique, the 
following model was used for the longitudinal motion’. 

u = (a^ - Aux ) - (q - Aq)w - gsin(0) (2.8) 

w = (a^ - Aa^) + (q - Aq)u + g cos(6) (2.9) 
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e = (q-Aq) 


( 2 . 10 ) 


where Aax , Aaz , Aq are the biases in accelerations along X and Z axis and pitch rate 
about Y axis, respectively. 

The observation equations used were^ 

ttjj, = K(^tan~*(Wjj/Un) + Aa 

where Aa is the bias in angle of attack (AO A) and Ka is the scale factor error in AO A. 
The following additional equations relating true airspeed at eg to true airspeed measured 
at the nose boom were used. 

Vn=A/K)^+(Wn)^ 

Un=u + (q-Aq)Zn 
Wn =w + (q-Aq)Xn 

where Xn and Zn are distances along X and Z axes from eg to pressure port. For the MLE 
program for the longitudinal case, the following set of unknown parameters was 
considered adequate^, 

=[Aax,Aaz,Aq,Aa,Ka] 

The initial conditions used for the aircraft states were the same as obtained from 
the measured time history. Further, it must be noted that since the three accelerometer 
locations were ofifeet from eg, suitable corrections were made to the measured 
accelerations. The accelerations a*, ay and a^ at c.g. are obtained from the accelerations 
axa, aya and a^ (measured at a point away from the eg ) through the following relations, 
ax = axa + (q^ +i^)Xa -(pq- f ) Y* -(pr+q )Za (2. 1 1 ) 
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ay = aya - (pq+r )Xa + (r^+q^)Ya -(rq- p )Za 


( 2 . 12 ) 


az = aza -(pr- q)Xa - (qr+ p ) Ya + (p^+q^)Za (2.13) 

where Xa , Ya and Za represent the offset distances of the accelerometers from the c.g. 
location. 

The onboard instrumentation scheme used for ICiran aircraft does not provide 
measurements of p, q and f. Therefore, the angular accelerations were obtained by 
numerical differentiation of the measured angular rates. Since the angular rates 
measurement were of reasonably good quality, it was expected that their differentiation 
will provide fairly accurate information about respective angular accelerations. 

For Lateral and Directional motion, the simplified model used w^as as follows^ 


u = -qw + (r - Ar)v 

(2.14) 

V = -(r - Ar)u + (p - Ap)w + (ay - Aay) 

(2.15) 

w = -(p - Ap) V + qu 

(2.16) 


where Ap, Ar, Aay are the biases in roll rate, yaw rate and lateral acceleration respectively. 
The observation equations used were 



where Un , Vn and Wn denote the velocity components along the three axes at the nose- 
boom, Ap is the bias in angle of sideslip (AOS) and Kp is the scale frictor error in AOS. 

For the lateral-directional motion, the following set of unknown parameters was 
considered adequate^. 
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e'^=[Aay,Ap,Ar,Ap,Kp] 

The MLE program was then run to estimate the respective correction factors for the 
longitudinal and lateral-directional motion. The cost function (refer Appendix 'B') was 
minimized by minimizing the error between the computed and measured observation 
variables. For longitudinal motion, the observation variables were true airspeed and 
angle of attack. For lateral-directional motion, the observation variables were true 
airspeed and angle of sideslip. 

The results of the data compatibility checks are as shown in Table 2.6a, 2.6b and 2.7 for 
the longitudinal and lateral-directional cases respectively. 


Table 2.6a Results of Data Compatibility Checks for FLT3 (Longitudinal Motion) 


SLNo- 

Parameters 

Units 


FLT3/TP13 

(Vi=180KIAS) 

1. 

Aax 


-1.0574 

-1.1398 



mu 

(0.0884)* 

(0.1853) 

2. 

Aax 






mmi 



3. 

Aq 

rad 


0.0429 





(0.0032) 

4. 

Aa 

rad 


0.0230 





(0.0052) 

5. 

K« 



0.6582 





(0.0195) 


* Cramer Rao Bounds 
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Table 2-6b Results of Data Compatibility Checks for FLT4 (Longitudinal Motion) 


SI. No. 

Parameter 

Units 

FLT4/TP2 

(Vi=160KlAS) 

FLT47rP8 

(V.=180KIAS) 

FLT4/rP15 

(Vi=220KIAS) 

1. 

Aax 

m/s^ 

-0.7117 

(0.1084)* 

-0.7732 

(0.0882) 

-0.8081 

(0.0623) 

2. 

Aax 

m/s^ 

-1.0578 

(0.3961) 

-1.0327 

(0.5422) 

3.3966 

(0.6192) 

3. 

Aq 

rad 

0.0367 

(0.0033) 

0.0292 

(0.0042) 

0.0076 

(0.0033) 

4. 

Aa 

rad 

0.0116 

(0.0029) 

0.0636 

(0.0040) 

0.0294 

(0.0045) 

5. 

Ka 


0.8238 

(0.0121) ! 

0.7019 

(0.0160) 

0.7544 

(0.0261) 


* Cramer Rao Bounds 


Table 2.7 Results of Data Compatibility Checks - (Lateral & Directional Motion) 


Parameter 

Units 

FLT5/TP1 
(V.= 160KIAS) 

FLT5/TP9 
(V.= 180KIAS) 

FLT5/TP17 
(V. = 200KIAS) 

FLT5/TP24 
(Vi = 220KIAS) 

Aay 

m/s^ 

-0.4332 

(1.1874)* 

-4.2584 

(1.9292) 

ne^Hi 


Ap 

rad/sec 


-1.2281 

(0.0248) 

-0.8129 

(0.0823) 

-0.6776 

(0.0560) 

Ar 

rad/sec 

IHH 

IHH 

IHHI 


AP 

rad 


-0.0569 

(0.0022) 


HilH 

Kp 




bIh 



* Cramer Rao Bounds 


Fig. 2.2a shows the match obtained between the measured and confuted angle of 
attack and true airspeed at the end of the data corrq>atibility check for the longitudinal 
motion. Similarly the match obtained between the measured and computed angle of 
sideslip and true ahspeed for lateral-directional motion are presented in Fig. 2.2b. 


These figures show that a good match is obtained between the conq)uted and 
measured observation variables. This implied that the correction fectors have been 
estimated with reasonable accuracy. From Table 2.6 and 2.7, it is observed that tte scale 
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Tas (m/s) 



Fig. 2.2b Data Compatibility Check - Comparision Of Measured And Computed 
Angle Of Sideslip(AOS) & True Airspeed (Tas) 
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factor for angle of attack (AO A) is of the order of 0.75 and that for angle of sideslip 
(AOS) is 0.85. This is also correlated physically by considering the location of 
AO A/ AOS vane. The vane is mounted in a plane just behind the blunt tip of the aircraft 
nose (the nose tip is blunt in shape because the landing light of the aircraft is fitted in the 
nose tip). The spindle on which the vane is moimted protrudes approximately 10 cm 
from the surface of aircraft body. This implied that the vane will be immersed in the near 
flow field of the aircraft nose. Hence, it is expected that the vane will e5q)erience flow 
interference as well as a certain amount of turbulent flow shed by the blunt nose tip, 
especially during maneuvers. The angle measured by the vane will be the 'local' angle. 
Therefore, the values of corrections factor Aa, AP, Ka and Kp , seems to be reasonable 
and acceptable. The correction fectors Aq and Ar were observed to be estimated with low 
values of Cramer Rao bounds and hence were acceptable. However, the values for 
correction frictor Aa* , Aay A^ and Ap are on the higher side then would normally be 
ejqpected. This could have been due to the feet that we have not estimated the initial 
conditions separately. Therefore, these correction factors were ignored. 

Thus, the raw flight data obtained from ASTE, lAF was corrected for the bias 
errors in an gular rates, the scale frictor and bias error for AOA and AOS, as well as for 
the ofifeet of accelerometers from the c.g. of the aircraft. The flight data obtained after 
applying these corrections is henceforth called as the 'corrected' flight data. Parameter 
estimation was then carried out on the corrected flight data as well as on the raw data 
using the Delta method. The details of the Delta method are now explained in the 
subsequent paragraph. 
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2.6 DELTA METHOD*’ 

The parameters occurring in the equations of motion of an aircraft represent 
partial derivatives of aerodynamic force and moment coefficients with respect to the 
corresponding motion or control input variables. The fundamental definition of partial 
derivatives is that they represent the change in the given function, when one variable is 
changed by a small amoimt, while the rest of the variables are held constant. Thus the 
parameters occurring in the equations of motion can be considered as variation in 
aerod)mamic coefficients due to a small variation in one of the motion or control variable 
(about its nominal value) while the rest of the variables are held constant at their nominal 
values. Let us now assume that the FFNN is already trained to m^ the network input 
variables to the given output variable. Figure 2.3a and 2.3b show a pictorial 
representation of the trained FFNN for longitudinal & lateral-directional cases. 


qc/2u 



INPUT VARIABLES OUTPUT VARIABLE 


Fig.23a Schematic Representation of the Trained FFNN For Longitudinal Motion 
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INPUT VARIABLES OUTPUT VARIABLE 

F^.2.3b Schematic Representation of the Trained FFNN For Lateral— Directional Motion 

Now only one of the network inputs is given a small (delta. A) perturbation at 
each time point, while ail others are held at their original values. Such a modified input 
file is then presented to the trained neural network to predict the perturbed value of 
aerodynamic coefficient at its output node. The difference in the predicted value of the 
aerodynamic coefificient from the value predicted for the original value of the input is due 
to the perturbation in the value of that particular network input variable. The difference, 
so calculated, in the value of the aerodynamic coefficient divided by the perturbation 
value yields the corresponding aircraft parameter. For example, say, a is perturbed by a 
small value Aa and the difforence observed in predicted value of Cm (for longitudinal 
case) fi)r a and a+Aa is ACm- Then ACjAa yields the stability derivative Cma- 
Similarly, the perturbation in fie and corresponding ACm observed wiU yield the control 
derivative Cmfie, where Cmfie=A CJAde. To avoid bias due one sided differences, 
motion and control variables are perturbed in both increasing and decreasing direction. 
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For example, if Cl values for a+Aa and a-Aa are predicted to be Cj^and Cl 
respectively, then 

CLa = (Cit-Ci;;)/2Aa 

It must be noted that for each of the perturbed value at each time point, there will be a 
corresponding predicted value of aerod5Tiamic coefficient and hence a different 
calculated value of the parameter estimate. Ideally speaking, all these values should have 
been identical so as to yield a single value of the parameter. However, in reality there is 
variation in the predicted values. The estimated values showed a near normal distribution 
and hence the mean value is used as the estimate and standard deviation about mean as 
the measure of accuracy of the estimates. 
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CHAPTERS 


RESULTS AND DISCUSSIONS 

3.1 SIMULATED FLIGHT DATA (LONGITUDINAL MOTION) 

The procedure as described in article §2.2 of Chapter 2 was adopted to generate 
the simulated flight data for longitudinal motion. The set of simulated flight data thus 
obtained will be referred to as the 'no noise' case, since they contain neither process nor 
measurement noise. This no noise data set was first used to estimate the longitudinal 
aerodynamic parameters (Cba » Ctq , Clsc , Cma , C^q and CmSe) by the Delta method. 

As required by the Delta method, the network input training file was prepared 

QC 

consisting of a, — and fie. The network output file consisted of the aerodynamic 
2u 

coefficient Cl or Cm . The FFNN was trained using these input-output file for different 
architectures of the FFNN. The training function used for training of the FFNN was the 
"Trainrp" function firom the Neural Network Toolbox of Matlab 5.3. The transfer 
function used for the neurons of the hidden layer of the FFNN was the standard "Log 
Sigmoidal" function from the Neural Network Toolbox of Matlab 5.3 (refer Eq. (A.1) of 
Appendix 'A"). The different architectures of the FFNN were obtained by varying the 
network parameters, such as learning rate, number of neurons in the hidden layer and 
logistic gain etc*'*. The architecture which gave the minimum "Mean Square Error (MSE)' 
was fi-ozen for the prediction phase*'* (refer Eq. (A.2) of Appendix 'A*). The table below 
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shows the range of values of the network parameters used for training the FFNN to map 


qc 

inputs (a, — ,oe) to output (Cl or Cm 
2u 


)• 


Network Parameter 

Range 

Logistic gain 

0.1 to 0.9 

Number of hidden layers 

1 

Number of iterations 

5000 

Learning rate 

0.1 to 0.9 

Number of neurons 

4 to 8 

Performance criteria 

Mean Square Error 


The number of iterations were kept fixed at 5000 because it was found that the 
decrease in MSB beyond 5000 iterations was marginal. 

Once the training phase is over, the Delta method was then used to estimate the 
aerodynamic parameters^^’^^’^^’^*. The values of the parameters estimated by applying the 
Delta method on the no noise simulated flight data (longitudinal) is presented in third 
column of Table 3.1. On conq)aring these estimates with the true values (used for 
generating the simulated data) given in column 2 of Table 3.1, it is seen that the estimated 
values are in good agreement with the true values. Further, the accuracy of the estimates 
are high as they have very low standard deviations. 

A study was also carried out to see the effect of measurement noise on the 
estimates of longitudinal parameters. Accordingly, simulated pseudo random noise of 
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varying intensity was added to the motion variables (a, -^) and the aerodynamic 

2u 

coefficients (CLjCm-)- The noise was simulated by generating successive uncorrelated 
pseudo random numbers having normal distribution with zero mean and assigned 
standard deviation, the standard deviation corresponds approximately to designated 
percentage (in the present work 5%, 10% and 15% were used) of the maximiun 
amplitude of the corresponding variable. The results of the parameter estimation for 5%, 
10% and 15% noise cases are presented in columns 4, 5 and 6 of Table 3.1. As seen from 
the Table 3.1, the Delta method yields good results even in the presence of noise. This 
implied that the addition of noise had a very marginal effect on the estimated parameters. 
In fact, except for Ctse and Cmq , the rest of the parameters have been estimated very well 
even in presence of 15% noise (which is unusually high level of noise). This only goes to 
reaffirm the robustness of the Delta method which is able to extract the parameters 
reasonably well even from noisy data**. In real life situations, the flight data will be 
noisy and the Delta method will be capable of handling it. 

In order to check the accuracy of the estimates, the aerodynamic coefficients 
(CljOii.) were estimated by feeding the values of the estimates obtained for the 15% 
noise case into the following equations. 

CL=Cu+CLaa + CLqr^l+CL5eSe (3.1) 

Cm = Cmo + Cma** I" 2u J (^-2) 

The values of the parameters Clo and Cmo were taken as zero in the above Eq. (3.1) and 
(3.2). Fig. 3.1 shows a conq)arison of the plot of the measured (simulated) and estimated 
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Cm 



Fig. 3.1 Plot of Measured and Ej 
F or Simulate 
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aerodynamic coefficients. From the Fig. 3.1, it is seen that the estimated coefficients lie 
very close to the corresponding measured coefficients. This implied that the estimates 
obtained from the Delta method are fairly accurate, even when estimated from flight data 
with unusually high measurement noise (15%). 

3.2 REAL FLIGHT DATA (LONGITUDINAL MOTION) 


3.2.1 General 

The real flight data for the longitudinal motion was generated from two flights, 
namely, FLT 3 and FLT 4, as discussed in article §2.4 of Chapter 2 (refer Table 2.5). 
From these two flights a total of eight data sets were obtained at different speeds. These 
are as follows : three data sets each for Vi = 160 KIAS and Vi = 180 BQAS and two data 
sets for speed Vi = 220 KIAS. The test input used for all the eight data sets was the 
3-2-1 -1 5e input. It must be noted that all the eight test inputs were not exactly identical. 
This is because firstly they were given by different test pilots and secondly four inputs 
began with positive elevator deflection and four with negative elevator deflection. Using 
the nomenclature system described in article §2.4 of Chapter 2, these eight real flight data 
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sets will hereafter be referred to as given below: 


Data Set 

Nomenclature 

1. 

FLT37tP 2 
(Vi =160 KIAS) 

2. 

FLT3/TP13 
(Vi = 180 KIAS) 

3. 

FLT4/TP1 
(V.= 160 KIAS) 

4. 

FLT4/TP2 
(V.= 160 KIAS) 

5. 

FLT4/TP7 
(V,= 180 KIAS) 

6. 

FLT4/TP8 
(V,= 180 KIAS) 

7 

FLT4/TP14 
(V, = 220 KIAS) 

8. 

FLT4/TP15 
(V. = 220 KIAS) 


It can be seen from the above, that the maneuvers have been repeated for a given 
speed. T his was done to ensure repeatability of data as well as to provide a sufficient 
database for analysis. 

The longitudinal parameters that are estimated in the present work are Cta , Ctq , 
CL8e , Qna , Qnq and C^se - Of thcsc, it is wcll loiown that Ctq and Clsc are weak 
derivatives, in the sense that, even with a large change in Ctq amd Ctse, there will be little 
or no appreciable change in Cl. In other words, the aerodynamic coefficient Cl is less 
sensitive to the longitudinal parameters Crq and Clsc- The major contribution to Cl is 

from Clk. 

3.2.2 Problem of Cl Obtained From az 

As was discussed in preceding article §3.1, for training of the FFNN, the fest step 

QC 

was to prepare the network input file consisting of a, and 5e variables. In order to 
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cater for the thrust effect, instead of a, q and 6e, the input variables a-a* , q-q* and 6e- 
5e have been used, where a , q and 6e are the trim values for the given test point. By 


doing so, the perturbed a, q and 6e have been modeled along with the perturbed Cl or 
Cm-. Accordingly, the output file was prepared with perturbed Cl or Cm- These 
coefficients were calculated as follows 


Cl « (3^2: V 

pu^S 


(3.3) 


CL=^[d-q] 

puS 


(3.4) 


C„,= [2qIy]/[pu2Sc] 


(3.5) 


Cl was calculated using as well as d from Eq. (3.3) and (3.4) respectively, a was 

obtained by the numerical differentiation of a. Therefore, a is boxind to contain errors 
due to numerical differentiation. The values of lift coefficients (Cl ) obtained from both 
the equations were compared and it was fovmd that the Cl calculated from showed very 
sharp and wild peaks. This was found to be so for all the test points analysed. In order to 
observe the relationship between Cl and a, the lift coefficient (Cl) obtained fi)r simulated 
data was plotted against a. For real flight data, both the lift coefficients (Cl from a^ and 

a) were plotted against a. These three plots are presented in Fig. 3.2. As seen from the 
Fig.3.2, the plot of Cl versus a for simulated data revealed that aU the points feO neatly 
mto a very narrow and linear band having a positive slope. However, the plot of Cl 
( from az) versus a, for real flight data, showed a large spread of the data points. In 

contrast, the plot of Cl (from a) versus a showed that although the points were spread in 


34 








a relatively broad band, yet clearly there was a linear pattern with positive slope, almost 
similar to that seen for simulated flight data. At this stage, it was conjectured that 
probably there was some problem in az recording. This was subsequently proved to be 
true after examining the values of the aerodynamic parameters estimated from Cl 
( calculated from az) using the Delta and MLE methods. The values of the estimates were 
found to be widely varying with very high standard deviations. Hence, it was decided not 
to use Cl calculated from az for estimation process. All the estimates for longitudinal 
parameters have been obtained using Cl calculated from d in the present work. 

3.2.3 Procedure Used For Estimating Parameters From Real Flight Data 

In the present work for the longitudinal motion, apart from estimation by the 
Delta method, the parameters were also estimated using MLE algorithm. This was done 
to provide a basis for comparison for the Delta method estimates with the MLE estimates. 
The MLE method requires a priori fix of the aerodynamic model. It is well known that 
the MLE method is sensitive to the initial guess values used. In absence of realistic guess 
values (based on analytical method or wind tunnel tests), initially arbitrary guess values 
were used as starting values. However, frequently it was observed that the MLE method 
with these guess values was not able to estimate large number of parameters to the 
desired accuracy. In order to run the MLE programs, the estimated values obtained from 
the Delta method of FFNN were used as the initial guess values for the MLE algorithm. 
It was observed that MLE estimates now were of desired accuracy in terms of sign and 
magnitude. 
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From Ref. 20, it was reported that FFNN are capable of extracting parameters 
from raw flight data containing bias and scale fector errors. In the present work, the 
parameters were estimated from raw flight data using the Delta and MLE method. 

It is well known that in order to improve the accuracy of the estimates^, data 
compatibility checks must be done on the raw flight data. Accordingly, data 
compatibility checks were applied (refer Chapter 2, article §2.5) to the raw data and 
'corrected' flight data was obtained. The parameters were then estimated from the 
corrected flight data using both the Delta and MLE methods. 

A systematic procedure, as given below, was followed to estimate the 
aerodynamic parameters from real flight data (longitudinal) for both the raw and 
corrected flight data using the Delta and MLE method. 

1. Aerodynamic parameters (Cta , CLq , Clsc , Cma , CnM, and Cmge ) are estimated by the 
Delta method xxsmg raw flight data. 

2. Using the above estimates as initial guess values, the MLE program is executed to 
estimate aerodynamic parameters (Clo , Ctq , Ctse , Cma > Cmq and Cmse ) from raw 
flight data. 

3. Next, the raw data is preprocessed using the correction factors given in Table 2.6 and 
2.7 in Chapter 2. In other words, data conqiatibility checks are done on the raw flight 
data. This flight data is referred to as 'corrected' flight data. 

4. Steps 1 and 2 are repeated for the corrected flight data. This is done to estimate the 
aerodynamic parameters using corrected data. 

The typical time history plots obtained for raw and corrected flight data following the 
above procedure are presented in Fig. 3.3 and 3.4. 
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3.2.4 Parameter Estimates From Raw Flight Data 

As explained above, the first step was to estimate the parameters from raw flight 
data using the Delta method. Accordingly, the aerodynamic modeling using different 
architectures of FFNN was carried out with network input training file having variables 
a-a* , q-q and 6e-5e* and output file having aerodynamic coefficients Cl (obtained from 
d ) or Cm . The architecture that gave the minimum MSB was frozen for the prediction 
phase. Fig. 3.5a and 3.5b show the plot of actual and predicted (after tr aining ) 
aerodynamic coefficients (Cl ,Cm.) obtained from the raw data. From the Fig. 3.5a and 
3.5b, it is seen that good training of the FFNN has taken place. In other words, the FFNN 
is able to map the input-output relationship accurately. This framed FFNN represented a 
reasonably good aerodynamic model for the aircraft under study. 

The Delta method was then used to estimate the parameters from raw flight data. 
The estimates for the longitudinal motion obtained from the Delta method of FFNN for 
the raw flight data are presented in Table 3.2a and 3.2b for flights FLT3 and FLT4 
respectively. 

As mentioned in step 2 above, for con^arison purpose, the parameters were also 
estimated by the MLB method from the raw flight data, using the Delta method estimates 
as initial guess values. The results of the MLB method are presented in Table 3.3a and 
3.3b for flights FLT3 and FLT4 respectively. 

From a study of results of the Delta method for raw flight data (Table 3.2a and 
3.2b), it is observed that all the strong derivatives (Clo , Cnw , Cmq and CniSc ) have been 
estimated with reasonably good accuracy. Further, the estimates from the Delta method 
(Table 3.2a and 3.2b) for all the test points at different speeds show a feirly consistent 
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trend and increase the level of confidence in the results. From the first row of Table 3.2a 
and 3.2b, it is also seen that the standard deviation obtained for Clo is slightly on the 
higher side for all the test points, ranging from 0.2674 for FLT4/TP14 to 0.5841 for 
FLT4/TP1 . Although the standard deviation for Clo are consistently on the higher side, it 
must be remembered that this parameter has been estimated from Cl obtained from d , 
with all its inherent errors due to numerical differentiation of a. Seen against this 
background, it can said that the value of Cta estimated using the Delta method are 
reasonably good, despite the slightly higher standard deviation. This reaffirms, as 
reported in Ref. 18, that the Delta method can be advantageously applied to estimate 
parameters from raw flight data, which may possibly contain scale factors and bias errors. 
This is a distinct advantage of the FFNNs. In fact, for future application one can quickly 
set up the FFNN and run it to have a quick look at the values of the parameters, if so be 
desired. Then subsequently, with little additional fine tuning of the network, reasonably 
good estimates can be obtained without the requirement of the preprocessing of the flight 
data. 

As seen from the Table 3.3a and 3.3b, the MLE estimates obtained for the strong 
derivatives are in close agreement with those obtained from the Delta method. It must be 
stressed that to initiate MLE algorithm, parameters estimated by the Delta method were 
used as initial guess values for the derivatives. However, as observed fix)m the Tables 3,2 
and 3.3, the weak derivatives Cl, and Clsc have not been estimated very well by either 
the Delta or the MLE method. 

After estimating the parameters from the raw data, the next step was to check the 
accuracy of the% estimates. This can be done by finding out the estimated total 
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aerodynamic coefficients (Cl, Cm) using the parameters estimated above and comparing it 
with the actual (measured) coefficients ( obtained from Eq. (3.4) and (3.5)). The 
estimated aerodynamic coefficients were calculated by feeding the parameters estimated 
above into the Eq. (3.1) and (3.2). However as seen from the Eq. (3.1) and (3.2), the 
parameters Clo and Cmo are not known, as they have not been estimated separately. 

In order to compute Clo and Cmo from the raw flight data, the following procedure 
was adopted**. Firstly, using the estimates from the Delta method, the estimated Cl and 
Cm are calculated from the Eq. (3.1) and (3.2) by ignoring the terms, Clo and Cmo - These 
estimated total aerodynamic coefficients (Cl, Cm) do not carry any information about 
parameters Clo and Cmo- However, the actual (measured) total aerodynamic coefficients 
of the aircraft carries the complete information, including Clo and Cmo (refer Eq. (3.4) 
and (3.5) given earlier). The above estimated coefficients were subtracted from the 
respective measured coefficients to obtain the values of Clo and Cmo**- However, since 
there are a number of time steps, the parameters Clo and Cmo have a value for each time 
step. The mean of this distribution of Clo and Cmo is taken to arrive at a single value. 
The results obtained for the different test points are shown m Table 3.4 . 

Figures 3.6, 3.7 and 3.8 show the comparative plots of the actual and estimated 
total aerodynamic coefficients (obtained from raw flight data) for three typical test points, 
after adding the parameters Clo and Cmo- These figures show that a ftiirly good match is 
obtained. This implied that the parameters have been estimated reasonably well at each 
test point fi)r different speeds. 

In order to further cort^are the values of the estimates obtained from the Delta 
and MLE methods from raw data, it was decided to compare the estimated and measured 
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responses. Accordingly, the set of parameters estimated using the Delta and MLE 
methods (from the raw flight data) were used to solve the longitudinal equations of 
motion (refer Eq. (2.1) and (2.2) of Chapter 2) . As stated earlier, these equations were 
solved using the fourth order Runga-Kutta method to get the longitudinal response of the 
aircraft in terms of a and q for each test point. 

The estimated responses were plotted along with the actual response of the 
aircraft. These plots are presented in Fig. 3.9, 3.10 and 3.11 corresponding to test points 
FLT3/TP2, FLT3/TP13 and FLT4/TP14 respectively. These plots reveal that a better 
match is obtained for FLT3/TP2 and FLT3/TP13 when using parameters estimated from 
the Delta method. This implied that the estimates from the Delta method are feirly 
accurate and that these values could be used to construct the mathematical model to 
generate the flight responses of the aircraft. However, in case of test point FLT4/TP14 
(Vi=220 KIAS) it is observed (Fig. 3.11) that the estimated response shows a slight time 
period shift as con:q)ared to the actxial response. One of the possible reasons could be the 
incorrect estimate of Clo> and Cnw obtained for this test point (refer column 4 of Table 
3.4). 

After completion of analysis of the raw flight da t a, the corrected flight data was 
taken up for parameter estimation (Step 3 of procedure explained at the beginning of this 
paragraph). It would be of interest to see how the results improve after ^plying the data 
compatibility checks to the raw flight data. 

3.2.5 Parameter Estimates From Corrected Flight Data 

The same procedure as described earlier for the raw flight data, was utilized here 
for the corrected fl ight data in order to obtain the parameter estimates from the Delta and 
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Fig. 3.9 Comparision of Actual And Estimated Response 
For FLT3fTP2,Vi=160 KIAS For Raw Data 
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Fig. 3.10 Comparision of Actual And Estimated Response 
For FLT3n'P13,Vi=180 KIAS For Raw Data 
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Fig. 3.1 1 Comparision of Actual And Estimated Response 
For FLT4/TP14,Vi=220 KIAS For Raw Data 




MLE method. Firstly, the aerodynamic modeling using the FFNN was carried out for 
the corrected flight data. Figures 3.12a and 3.12b show the plot of actual and predicted 
(after training) aerodynamic coefficients (Cl and Qn ) obtained from the corrected flight 
data. From Fig. 3.12, it is seen that the FFNNf has been well trained. The Delta method 
was then applied to estimate the parameters. The results of the Delta method for 
corrected flight data are presented in Table 3.5a and 3.5b for flights FLT3 and FLT4 
respectively. 

The results of the MLE method for corrected flight data are presented in Table 
3.6a and 3.6b for flights FLT3 and FLT4 respectively. From Table 3.5a and 3.5b, it is 
seen that the strong derivatives have been estimated reasonably well by the Delta method. 
The estimates for CLa , Qna , Qnq and CmSe (Table 3.5a and 3.5b) from the Delta method 
are fairly constant for all the test points at different speeds and show little variation. 
However, a few values do not fit into the pattern, for exan:q)le Cma has been obtained as - 
0.2463 for the data set of FLT4/TP8 (refer Table 3.5b), which is quite low. Further, it is 
observed that, with increase in speed, the value of CmSe decreases slightly, for example, 
CmSe = -0.7397 at V, = 160 KIAS and CmSe = -0.5526 at Vi = 220 KIAS ( refer last row of 
Table 3.5b). Strictly speaking, for the range of Mach number covered by different test 
points (from 0.3 M to 0.38 M) , the estimates slK>uld have been constant. The variation in 
the values of the estimates could be due to the feet that all the elevator inputs given 
during the maneuvers were not the ideal 3-2-1-1 5e ii^iut. 

On comparing the results of the Delta method for corrected flight data (refer 
Table 3.5a and 3.5b), with those for raw flight data (refer Table 3.2a and 3.2b), a few 
peculiarities were observed. The estimates of Cta » Cniq and CmSe » from the Delta n^thod 


54 



21 


Time (X 0.0625 second) 

Fig. 3. 1 2a Plot Of Actual And Predicted (After Training) CL For Test Point 
FLT3/TP2,Vi=160 KIAS For Corrected Data 



match feirly well between the corrected and raw flight data. However, the estimate for 
Cma appears to be suppressed for the corrected flight data (refer Table 3.5a and 3.5b) and 
consistently lower values of Cma are obtained here. 

As seen from the Table 3.6a and 3.6b, the MLE estimates from the corrected 
flight data are also fairly consistent and reliable. A comparison of the estimates from the 
Delta method (Table 3.5a and 3.5b) and the MLE method (Table 3.6a and 3.6b) shows 
that the estimates from both the methods are in close agreement. The Delta method can 
therefore be successfully employed for parameter estimation from the corrected flight 
data. 

The accuracy of the parameters estimated from the corrected flight data by the 
Delta method was checked in a similar fashion as was done for the raw flight data. The 
estimates for the parameters Clo and Cmo obtained for the corrected flight data are 
presented in Table 3.7. 

The estimated aerodynamic coefficients (Cl , Cm ) were then found out, using the 
estimates for corrected flight data obtained by the Delta method. Figures 3.13, 3.14 and 
3.15 show the plots of the actual and estimated aerodynamic coefficients (Cl and Cm) for 
corrected flight data for three typical test points. These plots show that a reasonable 
match is obtained between the measured and estimated aerodynamic coefficients. 

The accuracy of the estimates from the corrected flight data using the Delta and 
MLE methods were also checked by con:q)aring the measured and estimated reqionses. 
The method used to find out the estimated responses was similar to that followed for the 
raw fli ght data. The estimated responses for the corrected flight data were then plotted 
along with the measured response of the aircraft as shown in Fig. 3.16, 3.17 and 3.18. 
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Fig. 3. 1 6 Comparision of Actual And Estimated Response For 
FLT3/TP2,Vi=160 KIAS For Corrected Data 
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Fig. 3. 1 7 Comparision of Actual And Estimated Response For 
FLT3/TP13,Vi=180 KIAS For Corrected Data 









Fig. 3. 1 8 Comparision of Actual And Estimated Response For 
FLT4n‘P14, Vi=220 KIAS For Corrected Data 
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From these figures, it is observed that, in general, for the corrected flight data, the match 
between the estimated and measured response are not as good as that for raw flight data 
(Fig. 3.9, 3.10 and 3.11). However, in case of the corrected flight data, the estimated 
response using the estimates obtained from the Delta method hes very close to the 
measured response. Further, at peaks it is observed that there are small differences in 
absolute values of a and q (Fig. 3.16, 3,17 and 3.18). Also, the time period for the 
estimated response shows a slight increase (implying a decrease in short period 
frequency) as compared to the measured response. 

On completion of the parameter estimation process, the distribution of the 
parameters estimated by the Delta method from the raw flight data was plotted and is as 
shown in Fig. 3.19. As seen from Fig. 3.19, the estimates obtained from the raw flight 
data using the Delta method are fairly accurate and reliable. This is further corroborated 
by the fact that, when these Delta method estimates are used as initial guess values in 
MLE method, the MLE method also gives feirly accurate and consistent results. 
However, the estimates obtained from the corrected flight data using the Delta method 
can be improved further by carrying out a rigorous data conq^tibility check (using the 
complete six degree of freedom equations of motions®). In the present work, this could 
not be done because the pitch angle (0) transducer was unserviceable. Hence a reliable 
recording of pitch angle (0) was not available. Therefore, a simplified set of longitudmal 
equations of motions (refer Eq. (2.8), (2.9) and (2.10) given in article §2.5 of Chapter 2) 
were used for data com patib ility checks. Further, it was assumed that the mitial 
conditions for the data compatibility checks were as recorded during the beginning of the 
maneuver. Probably, the estimation of the initial conditions separately may have helped 
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Fig. 3. 1 9 Plot Showing The Distribution Of The Parameters Estimated Using 
The Delta Method From The Raw Flight Data 
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in improving the accuracy of the correction factors. And lastly, one of the observable 
variables used was true airspeed, which in turn was calculated from the dynamic pressure 
recorded. The transducer used being a pressure transducer is bound to have lags/time 
delays in its recording. In spite of these discrepancies, it is felt that the match obtained 
between the estimated response obtained from the Delta method and measured response 
for the corrected flight data is reasonably good. 

3.3 REAL FLIGHT DATA (LATERAL-DIRECTIONAL MOTION) 

3.3.1 General 

The real flight data for the lateral-directional motion was generated from the last 
flight of the test program, i.e. from FLT5, as discussed in article §2.4 of Chapter 2. 
Based on Ref. 9 and 19, it was decided to excite the lateral-directional response of the 
aircraft by giving separate 8a and 5r inputs in different test points. This meant that ft)r 
the first test point when 8a input is being given, 8r is held constant at its trim position and 
vice versa for the second test point. The dutch roll in Kiran aircraft is easily excited by 
giving full 8r doublet input alone. Hence for 8r iiq)uts, it was decided that only the 
doublet form of input will be used. However, for aileron (8a) inputs, the following input 
form were tried out; 3-2-1-1 H^)ut of time period of 7 seconds, 8a doublet and roll 
maneuvers’. Further, the flight data for lateral-directional motion was generated at four 
different speed, namely, Vi = 160, 180, 200 and 220 KIAS. 

From aerodynamics, it is well known that the rolling moment coefficient (Cj), 
yawing moment coefficient (Cn) and side force coefficient (Cy) are a function of F(p, r, p, 
8a, 8r). Hence in this case, the number of parameters to be estimated (in order to 
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determine the function 'F'), becomes large and any parameter estimation process becomes 
an involved exercise. In the present work, we have atten^ted to estimate the following 
parameters : Qp , Cir , Cip , Cisa , Ci8r 5 Cnp 5 Cnr 5 CnP ? CnSa » CnSr 5 f-Vp ’ ’ CySa •> 

CySr • 

3.3.2 Problem Faced During Parameter Estimation Using MLE Method 

The success in estimating the longitudinal parameters using the Delta method 
embolden us to work directly with the real flight data for this case. Accordingly, no work 
was done on simulated data for lateral-directional case. Further, like in the longitudinal 
case, after completion of estimation from the Delta method, it was also attempted to 
estimate the lateral-directional parameters from the real flight data using the MLE 
method. Here, again, the Delta method estimates were used as the initial guess values to 
initiate the MLE program. But, despite rigorous efforts, we could not succeed with the 
MLE method within the short span of time available to us. The reason for this was very 
simple. In the MLE algorithm, all the 15 parameters were being estimated 
simultaneously. Hence, it was found that the MLE method had become very sensitive to 
the initial guess values used. It was observed that for any given set of initial guess values 
used, the MLE method correctly estimated only a few parameters. A subsequent run, 
with a minor change in some of guess values, resulted m few other parameters being 
estimated correctly. Thus convergence of the MLE method was found to be a direct' 
function of the initial guess values. After considerable discussions, it was suspected that 
maybe the MLE routine written in house was incorrect To prove that this was not so, a 
set of simulated data was generated for the lateral-directional motion using the Delta 
method estimates (as true values of the aerodynamic parameters) and 8r doublet input as 
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given in actual flight. The Eq. (2.3) to (2.7) as given in Chapter 2 were used fijr 
generating the simulated data. Now, the same MLE routine was run using the Delta 
method estimates as initial guess values for the parameters. The results of this small 
exercise is placed in Table 3.8. It must be noted that the Delta method estimates used as 
initial values in the MLE routine were suitably altered by a small random number (refer 
column 3 of Table 3.8). It was found that the MLE routine quickly converged in the sixth 
iteration itself and gave us all the correct estimates with minimum cost function being 
foxmd as 2.1456x10’*'*. At this stage, it was even more keenly appreciated how an 
alternate method, which does not require initial guess values to start the estimation 
process, would be extremely helpful and useful for such similar situations in industry. 
No further attempts were made with the MLE method for real flight data and we proceed 
to discuss the Delta method results. 

3.3.3 Results of Parameter Estimates From Flight Data Generated by Individual 
5a and 5r Inputs 

The Delta method was used to estimate parameters from the raw flight data as 
well as from the corrected flight data as was done for the longitudinal case. Each test 
point was analyzed separately. Accordingly, for the 5r doublet input test point, the 

network training files containing — , — , p and 5r as input variables and Ci , Cn or Cy as 

2u 2u 

the output variable were prepared. Similarly, for the 5a iiqjut test point, the network 

input file contained — , — , p, and 5a with output file having Ci , Cn or Cy . The actual 
2u 2u 

(measured) aerodynamic coefficients Cj and Cn were calculated using Eq. (2.4) and (2.5) 


68 


as listed in Chapter 2. The actual aerodynamic force coefficient Cy was calculated from 
the lateral acceleration ay as follows: 

a = -^Cy + (pq + f)Xa - (r^ + + (rq - p)Za (3-6) 

mg ^ 

Using such input-output files, the aerodynamic modeling using different 
architectures of FFNN was carried out. As explained earlier for the longitudinal case, the 
architecture of FFNN which gave the minimum training MSB was frozen for the 
prediction phase. The Delta method was then appUed to estimate the parameters. 

It was observed that, since the maneuvers were executed by giving separate 6r, 5a 
inputs, the estimates of the corresponding cross derivatives (Ct , Cisr » Cnp and CnSa ) 
obtained by the Delta method (for raw flight data) were showing lot of scatter and 
inconsistency. Also the standard deviation obtained in the above cases were very high. It 
was conjectured that this was due to the feet, that for Kiran aircraft, the individual control 
inputs (5a and 5r input separately) are probably not sufficient to the excite the complete 
lateral-directional motion of the aircraft. 

3.3.4 Revised Approach For Parameter Estimation 

From the above discussions, it is clear that the approach used so fer for the 
estimation process for lateral-directional motion needed to be revised. Several 
researchers had recommended that although the 5a and 5r input should not be given 
simultaneously, but giving the 2 inputs delayed from each other by a small time interval 
gave satisfactory results. However, the data available with us was for each input 5a and 
5r being given separately. After considerable debate, we decided to try out another 
approach, wherein we used the feet that for the Delta method, we can always concatenate 
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data from different time segments. It must be noted that this is not possible for the MLE 
method. In other words, the Delta method is independent of the time at which the data 
was generated. In this new approach; we concatenated the data from the two different 5a 
and 5r input test inputs for the same speed. It is important to remember here that the 
aerodynamic coefficients Ci , Q, and Cy must be individually calculated for each test 
point using all the relevant conditions for the given test points, such as weight, dynamic 
pressure etc. The steps involved using this new approach for parameter estimation by the 
Delta method were as follows: 

1) . For each of the test points for 5a and 5r inputs (which were given separately) for the 

same speeds, the separate files containing input variables — , — , P, and 5r and 5a with 

2u 2u 

corresponding output variable Q , Cn and Cy were prepared. 

2) Now the two files were merged by putting the 5a input file at the end of the 5r input 
file, to form a combined training file. Such a combined file would carry the complete 
lateral-directional response of the aircraft being excited by 6r and 5a inputs. 

The next issue to be resolved was, as to which of the 5a inputs should be used to 
concatenate the data. The real flight data was generated for three different types of 5a 
input. A brief study revealed here that the 6a doublet input gave better estimates for 
rolling moment derivatives, when this test point was analyzed separately. 

Accordingly, it was decided that for the raw and corrected flight data, the aircrafr 
response obtained for separate 6r doublet input and 6a doublet input would be combined 
for each speeds. Thus four combined training files were prepared, one each for V,- = 160, 
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180, 200 and 220 BQAS. Accordingly, these four combined data sets have been 
designated as follows: 


Data Set 

Nomenclature 

1. 

FLT5/TP1 5 
(Vi= 160KIAS) 

2. 

FLT5/rP9 14 
(V,= 180KIAS) 

3. 

FLT5/TP17 22 
(V. = 200 KIAS) 

4. 

FLT5/TP24 28 
(V.== 220 KIAS) 


Figures 3.20 and 3.21 shows the typical time history plots of one such combined file for 
data set of FLT5/TP1_5 for raw and corrected flight data respectively. 

3.3.5 Parameter Estimates From Raw Flight Data Using Revised Approach 

Using the combined input-output training files, the aerodynamic modeling using 
different architectures of the FFNN for lateral-directional motion was carried out. The 
architecture of FFNN, which gave the minimum MSE, was frozen for the prediction 
phase. Figures 3.22a, 3.22b and 3.22c shows the plot of actual and predicted (after 
training) aerodynamic coefficients for a typical combined test point for raw flight data . 
A close match is obtained implying that the FFNN was now well trained. The Delta 
method is then used to estimate the parameters from the raw flight data. 

The results of the estimates from the Delta method for raw flight data using these 
four combined input files are presented in Table 3.9. From Table 3.9, it is seen that all 
the strong derivatives have been estimated very well with low value of standard 
deviation. The estimates are also observed to be feirly constant for all the test points at 
different speed. As seen from the first row of Table 3.9, the estimate of the roll damping 
derivative, Qp , has been obtained nearly constant as -1.1 but with slight higher standard 
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beta (radian) 



Time (X 0.0625 second) 

Fig. 3.20 Plot Of Time History For Combined Test Point 
FLT5/TP1-5, Vi=160 KIAS For Raw Data 
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Fig. 3.20 Continued 
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Fig. 3.20 Concluded 
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Fig. 3.21 Plot Of Time History For Combined Test Point 
FLT5n‘P1-5,Vi=160 KIAS For Corrected Data 
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Fig. 3.21 Continued 
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deviation, especially for test point FLT5/TP9_14 (Vj = 180 KIAS) \^*ere the standard 
deviation for Qp is 0.28. However, from the other three test points, the standard 
deviation is seen to be of the order of 0.18. This is considered acceptable, as it implied 
that the value of Cip hes between -1.298 to -0.902. 

Similarly, the strong control derivatives Cisa and CnSr have been obtained as -0.35 
and -0.12 with standard deviation of the order of 0.08 and 0.04 respectively (refer rows 
4 and 10 of Table 3.9). The dutch roll damping derivatives Cnr shows a slight scatter in 
its values for different speeds, from -0.92 at Vj = 160 KIAS to -0.55 at V, = 220 KIAS, 
but has been estimated in each case with a low standard deviation (refer row 7 of 
Table 3.9). The remaining strong derivatives Qp , Qp and Qp have been estimated 
consistently as -0.12, 0.22, -0.9 for all the test points with low standard deviation (refer 
rows 3, 8 and 13 of Table 3.9). 

The cross coupling derivative Qr estimated shows a slight variation, especially for 
the test point, FLT5/TP9-14 (Vj = 180 KIAS) with a correspondingly higher standard 
deviation (refer row 2 of Table 3.9). But the estimates for Q from the other three test 
points show good consistency with low standard deviation and lead us to place more 
confidence in these values. 

In feet, even the weak derivatives such as Qsr , Qp and Qsa have been estimated 
with remarkable consistency for different test points (refer rows 5, 6 and 9 of Table 3.9). 
Also, the derivative Qpr, except for test point FLT5/TP9_14, (Vj = 180 BOAS) has been 
estimated reasonably well (refer last row of Table 3.9). However the weak derivatives of 
sideforce coefficients, Qp , Qr and Qaa are found to be varying and inconsistent for each 
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of the test point analyzed and they have not been estimated very well (refer rows 11,12 
and 14 of Table 3.9). 

Overall it can be said that, except for these few weak derivatives, the results of the 
Delta method, as presented in Table 3.9, show good consistency among the estimates for 
different speeds. 

To check the accuracy of the estimates, the estimated total aerodynamic 
coefficients (Ci , €„ and Cy ) were compared with the actual (measured) coefiScients 
(calculated earlier from Eq. (2.4), (2.5) and (3.6)). The estimated aerodynamic 
coefficients were obtained by feeding the above estimates into the following equations. 

c 1 = Cip(^) + C|,(^) + CipP + CisaSa + Ci8.5r (3.7) 

~ ^np ^ ^ ^n5a^ (3-8) 

Cy = Cyp C";^) ^yr(^) ^ypP Cy5r^^ (3-9) 

These estimated toted aerodynamic coefficients do not carry any information 
about the parameters Cio , Cno and Cyo , since these paranaeters were not estimated 
separately. The same procedure as given in previous article §3.2, was used to estimate 
Cjo , Cno , Cyo. The estimates for the parameters Cio , Cno and Cyo obtained for raw flight 
data is presented in Table 3.10. 

Figures 3.23, 3.24, 3.25 and 3.26 show the conparative plots of the actual 
(measured) and estimated Ci , Cn and Cy for the raw flight data for all the four combined 
test points. From these figures, we see that a feirly good match is obtained for Cn and Cy 
for all the test points. However from Fig. 3.23, 3.24 and 3.25 (for Vj = 160, 180 and 200 
KIAS combined test points) for Ci , a small peculiarity is observed. It is seen that 
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although a good match is obtained between the measured and estimated Q over a large 
portion the figure, however, in a small time segment of the figure, the match is rather 
poor. This time segment corresponds to the beginning of the 8a doublet input. It can 
only be conjectured that this may have been due to the fact that 8a and 8r inputs were 
given separately in different test points (at same speed). In future, for lateral-directional 
motion, the 8r input should be given first and following 1-2 second delay, the 8a input 
should be given in the same test points. 

From the reasonably good match obtained between the actual and estimated 
aerodynamic coefficients for all the other test points (refer Fig. 3.23, 3.24, 3.25 and 3.26), 
it can be inferred that the set of parameters estimated by the Delta method from the raw 
flight data are fairly reliable and accurate. 

33.6 Parameter Estimates From Corrected Flight Data Using Revised Approach 

After conq)ietion of parameter estimation process using the raw flight data, the 
corrected flight data obtained after the data compatibility checks was taken up fijr 
analysis. During the data compatibility checks on the raw flight data, the correction 
factors given in Table 2.7 of Ch^ter 2 were used. The same procedure as was used for 
the raw flight data, was utilized here for parameter estimation by the Delta method. In 
this case also, the combined training files for the four speeds, Vj = 160, 180, 200 and 220 
KIAS were prepared. The actual aerodynanuc coefficient Q, Q and Cy were calculated 
using Eq. (2.4), (2.5) and (3.6). Next, the aerodynamic modeling using the FFNN was 
carried out for the corrected flight data. Figures 3.27a, 3.27b and 3.27c shows the plot of 
actual and predicted (after training) aerodynamic coefficients obtained for a typical 
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0.03 



Fig. 3.27a Plot Of Actual And Predicted (After Training) Cl For Combined 
Test Point FLT5/TP17-22,Vi=200 KIAS For Corrected Data 
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Time (X. 0.0625 second) 


Fig. 3.27b Plot Of Actual And Predicted (After Training) Cn For Combined 
Test Point FLT5/TP17-22,Vi=200 KIAS For Corrected Data 






Fig. 3.27c Plot Of Actual And Predicted (After Training) Cy For Combined 
Test Point FLT5n'P17-22,Vi=200 KIAS For Corrected Data 
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combined test point. From Fig. 3.27a, 3.27b and 3.27c, it is seen that good training of the 
FFNN has taken place. The Delta method was then used to estimate the parameters from 
this trained FFNN. 

The results of the estimation process using the Delta method for the corrected 
flight data is presented in Table 3.11. A study of these results shows that the strong 
derivatives, Cip , Qp and Qsa have been estimated reasonably well (refer rows 1 , 3 and 4 
of Table 3.11). For the control derivative CnSr , except for test point FLT5/TP24_28, 
(V i = 220 KIAS), a good match is obtained for the rema ining three test points (refer row 
10 of Table 3.11). In addition, the cross derivatives, Qpr and CnSa are estimated 
consistently with reasonable accuracy as indicated by the low standard deviation (refer 
row 5 and 9 of Table 3.1 1). In this case too, it is found that the estimate for Qr for test 
point FLT5/TP9-14, (Vj = 180 KIAS) is quite low with a high standard deviation (refer 
row 2 of Table 3.11). But the estimate of Ct for remaining three data sets matches very 
well, with low standard deviation being obtained in each case. 

The estimates for Qp , Cnp and Cyp derivatives shows good match for aU the test 
points analyzed for the corrected flight data. However, from Table 3.1 1, it is seen that for 
corrected flight data, the weak derivatives Cyp , Cyr and CySa have not been successfully 
estimated. 

On comparing the results, for the raw flight data (Table 3.9) with the corrected 
flight data (Table 3.11), it is observed that the estimates for derivative, C|p , Qp , Cipa , 
Cisr, Cnp, Qir, Cnp, Cn5a and C„8r are feirly consistent. However, for Cyp it is seen that 
slightly lower values are obtained for the corrected flight data. Thus, the Delta method is 
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capable of extracting parameters from the raw as well as from the corrected flight data for 
lateral-directional motion. 

In order to check the accuracy of the estimates obtained from the corrected flight 
data, the estimated and actual (measured) aerod 5 mamic coefficient,- Ci , Q, and Cy were 
compared. The parameters Qo , Cno and Cyo for the corrected flight data were estimated 
as was done for raw flight data. The results of the estimates for the parameters Cio , Cno 
and Cyo are presented in Table 3.12. 

Figures 3.28, 3.29, 3.30 and 3.31 show the plots of the actual (measured) and 
estimated lateral-directional aerodynamic coefficients for the corrected flight data for all 
the combined test points. These figures reveal that although a good match is obtained for 
Cn and Cy , in this case too, the same peculiarity for Cj is observed, as was seen for the 
raw flight data. In view of the good match obtained (refer Fig. 3.28, 3.29, 3.30 and 
3.31), it can be stated that the parameter estimated using the Delta method are feirly 
accurate. 

After completion of the parameter estimation process for the lateral directional 
case, the distribution of the parameter estimated using the Delta method from raw flight 
data was plotted and is as shown in Fig. 3.32. From the figure, it is seen that all the 
strong derivatives have been estimated with a very small scatter. 
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Table 3.1 Parameter Estimates From Simulated Data For Longitudinal Motion 


Parameter 
(Per Radian) 

True 

Value 

Estimates for V, = 160 KIAS 

NO NOISE 

5% NOISE 

10% NOISE 

15% NOISE 

Clo 

49 

4 8841 ■ 

(0 0476)* 

4 8239 
(0 0914) 

4 9258 
(0.0291) 


CLq 

10 0 

9 7981 
(0 0053) 

8.7800 

(0.0170) 

iSEmaliH 



04 

0 4136 
(0 0225) 

0.3808 

(0.0154) 



Cma 

-0 60 

-0 5991 
(0 0032) 

-0 6090 
(0 0075) 



Cmq 

-13.0 

-12.7540 

(0.0012) 

-12.7868 

(0.0024) 



Cni5e 

-0 68 

-0.6753 
(0 0044) 

-0.6849 

(0.0063) 

-0 6780 
(0.0182) 

-0 6640 
(0 0140) 


* Standard Deviation 
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Table 3.2a Parameter Estimates from Raw Flight Data Using Delta Method of 
FFNN for FLT3 


Parameter 

FLT3/TP2 

FLT3/TP13 

(Per Radian) 

(V.= 160KIAS) 

(V,= 180KIAS) 

Cha 

3 8822 

4.0894 


(0.5481)* 

(0.4596) 

CLq 

46 8668 

15 67932 


(0 3942) 

(0 2836) 

ClSs 

-0 6139 

-0 8556 


(0 9352) 

(0 4080) 

Qma 

-0.4689 

-04178 


(0 1410) 

(0.0867) 

Cmq 


-15 5018 



(0 0270) 

Cm5e 

-0.8649 

-0 8028 


(0.1676) 

(0 0514) 


* Standard Deviation 


Table 3.2b Parameter Estimates from Raw Flight Data Using Ddta Method of 
FFNN for FLT4 


Parameters 
(Per Radian) 

Estimated Parameters 


Vi = 16C 

)KIAS 

Vi =180 KIAS 

Vi = 220 KIAS 


FLT4/TP1 

FLT4/TP2 

FLT4/TP7 

FLT4/TP8 

FLT4/TP14 

FLT4/TP15 

Cloc 




4.2679 

(0.5045) 

4.2349 

(0.2674) 

4.1688 

(0.3853) 

CLq 


8.4112 

(0.6954) 


21.3603 

(0.1537) 

33.4108 

(0.3254) 

-164.9325 

(0.2874) 

ClSb 

ggH 

0.1017 

(0.4208) 


^bI 



Cma 


-0.4146 

(0.0370) 

ibh 

BiMI 

-0.3948 

(0.1089) 

-0.3984 

(0.0905) 

Ctnq 

-11.8338 

(0.0308) 

-12.9635 

(0.0105) 

-9.9484 

(0.0192) 


-8.8664 

(0.0554) 


Cm6e 

-0.6659 
(0 0623) 

-0.7258 

(0.0601) 

-0.5907 

(0.0958) 

-0.6738 

(0.0523) 

-0.5392 

(0.1778) 

-0.5556 

(0.0321) 


* Standard Deviation 
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Table 3.3a Parameter Estimates from Raw Flight Data Using MLE Method for 
FLT3 


Parameter 

FLT3/TP2 

FLT3/TP13 

(Per Radian) 

(V. = 160 KIAS) 

(V.= 180 KIAS) 

CLa 

4 1610 

4 6630 


(0 2978)* 

(0 5932) 

Cl^ 

117 3570 

241 2067 


(20 0271) 

(23 5703) 

ClSc 

2 1487 

4 6854 


(0 5618) 

(0 7942) 

Cma 

-0 5469 

-0.5766 


(0 01936) 

(0 0364) 

Cmq 

-12 0034 

-12.5565 


(1 2316) 

(1 9991) 

Cm5e 

-0 8063 

-0 7519 


(0 03669) 

(0.0399) 


* Cramer Rao Bounds 


Table 3.3b Parameter Estimates from Raw Flight Data Using MLE Method for 
FLT4 


Parameters 
(Per Radian) 

Estimated Parameters 


V.= 160 KIAS 

Vi = 180 KIAS 

V, = 220 KIAS 



FLT4/TP2 



FLT4/TP14 

FLT4/TP15 

Cha 


3 3498 
(0.2736) 


4.4009 

(0.4005) 

4.3089 

(0.5734) 

4.6888 

(0.1622) 

CLq 

92.0340 

(35.1748) 

61.2018 

(16.7365) 

8.1656 

(30.1511) 

-0.7120 

(38.6795) 

-33.0609 

(43.1463) 

5.3840 

(15.2477) 

ClSb 

|EQsnn|| 

1 1173 
(0.3621) 

-1.5487 

(0.7018) 

-0.1264 

(0.7371) 

-2.5254 

(0.9482) 


Cma 

-0.4786 

(0.0303) 

-0.5106 
(0 0157) 

Mil 


-0.4693 

(0.0275) 


Cmq 

-12 9971 
(1.8550) 

^^9 

nisi 

-16.1278 

(1.8551) 

-8.2094 

(2.0564) 

-17.6743 

(0.7868) 

Cm6e 

-0 6718 
(0.0362) 


-0.6879 

(0.0291) 

— 

-0.5764 

(0.0410) 

-0.6888 

(0.0168) 


*Cramer Rao Bounds 
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Table 3.4 Estimates of Parameters Cl© and Cm© from Raw Flight Data 

(Longitudinal Motion) 


Parameter 

Test Points I 


FLT3/TP2 

FLT3/TP13 

FLT4/TP14 


(V, = 160 KIAS) 

(V,= 180 KIAS) 

(V. = 220 KIAS) 

Clo 

0.1602 

0.1178 

0.0255 

Cmo 

0.0067 

0.0263 

0.0152 



Table 3.5a Parameter Estimates from Corrected Flight Data Using Delta Method of 
FFNN for FLT3 


Parameter 

FLT3/TP2 

FLT3/TP13 

(Per Radian) 

(V,= 160KIAS) 

(V, = 180 KIAS) 

CLa 

4 0101 

3.9773 


(0 6019)* 

(0 2218) 

CLq 

-119 0009 

-245 8390 


(0 3928) 

(0 2731) 

ClSb 

-2 2421 

-2.7677 


(0 7699) 

(0 7701) 

Cma 

-0 3684 

-0 3693 


(0.0451) 

(0.1226) 

Cinq 

-14 1881 

-15.9615 


(0.0125) 

(0 0201) 

Cni5e 

-0.7843 

-0.7824 


(0 0836) 

(0.0767) 


* Standard Deviation 


Table 3.5b Parameter Estimates from Corrected Flight Data Using Delta Method of 
FFNN for FLT4 


Parameters 

Estimated Parameters 


Vi = 160 KIAS 

Vi= 18 

OKIAS 

V = 220 KIAS 


FLT4/TP1 

FLT4/TP2 

FLT4/TP7 

FLT4/TP8 

ELT4/TP14 

FLT4/TP15 

CLa 

3 9771 
(0.5696)* 




4 0128 
(0.5414) 


Cu, 

-91.6920 

(0.2968) 

-121.348 

(0.3769) 

MM 

-150.211 
(0 2489) 

-126.1590 

(0.4613) 

-120.447 

(0.1464) 

ClSb 

BQg|l[|| 

-0.6605 

(1.4762) 

-2.0517 

(0.6930) 

Wmm 

-2.1410 

(0.2065) 

-0.9724 

(0.3360) 

Cma 

-0 3524 
(0 0975) 

-0.3447 

(0.0278) 

-0.3023 

(0.0585) 



-0.3410 

(0.0898) 

C„.q 


-15.2054 

(0.0053) 



mi 

-12.4383 

(0.0440) 

Cm6e 


MM 

-0.5797 

(0.0528) 

-0.6133 
(0 0977) 

-0.5754 

(0.0511) 

-0.5526 

(0.0653) 


* Standard Deviation 
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Table 3.6a Parameter Estimates from Corrected Flight Data Using MLE Method 
forFLTS 


Parameter 

FLT3/TP2 

FLT3/TP13 

(Per Radian) 

(V, = 160 KIAS) 

(V,= 180 KIAS) 

CLa 

4.1527 

4 3309 


(0 4227)* 

(0 3060) 

Clq 

-25 6215 

4 3663 


(34 2281) 

(28 9447) 

ClSb 

2 5199 

49117 


(1 0732) 

(0 7745) 

Cma 

-0 4012 

-0.3379 


(0 0171) 

(0 0119) 

Cmq 

-12.2771 

-14.5814 


(1 5311) 

(1 2183) 

Cm6e 

-0 8202 

-0 7771 


(0 0393) 

(0 0262) 


* Cramer Rao Bounds 


Table 3.6b Parameter Estimates from Corrected Flight Data Using MLE Method 
for FLT4 


Parameters 

Estimated Parameters 


Vi = 160 KIAS 

Vi = 180 KIAS 

Vi = 220 KIAS 


FLT4/TP1 

FLT4/TP2 

FLT4/TP7 

FLT4/TP8 . 

FLT4/TP14 

FLT4/TP15 

CLa 

2.5457 

3.3539 

4.0731 

4.0638 

4.5975 

5.2662 





(0.7153) 

(0.5025) 

CLq 

18 9475 

-41.3614 

-136 032 

23.3044 




(21 5061) 

(39.9195) 

(41 2101) 




-0.7376 

0.9846 

-2.2189 


-4.2933 

-2.6445 

(0 7326) 

(0.5009) 

immm 


(1.6929) 

(1.4310) 

Ciuot 

-0 4061 

-0 4181 

-0.3155 

-0 3909 

-0.3518 

-0.2615 


(0.0259) 

(0.0128) 

(0.0146) 

(0.0236) 

(0.0243) 


Cmq 

-12 0530 

-11.1397 

-12.3961 

-13.5502 

-10.4634 

-17.2019 

(1.7880) 

(0.9603) 

(1.3975) 

(1 9995) 


(1.8072) 

Cm5e 

-0.6263 


-0.6876 

-0.6472 

-0.6100 


(0.0392) 

(0.0251) 

(0.0299) 

(0.0420) 

wmsm 

(0.0285) 


*Cramer Rao Bounds 
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Table 3.7 Estimates of Clo and Cm© from Corrected Flight Data (Longitudinal 
Motion) 


Parameter 

Test Points 


FLT3/TP2 

FLT3/TP13 

FLT4/TP14 


(V,= 160 KIAS) 

(V,= 180 KIAS) 

(Vi = 220 KIAS) 

Clo 

0.1192 

0.0751 

-0.0064 

Cmo 

0.0023 

0.0078 

0.0172 



Table 3.8 Parameter Estimates From Simulated Data (Lateral-directional Motion) 
Using MLE Method 


Parameter 

True Values of 
Parameters 

Starting Guess Values 

Estimates 

Cip 

-1 2032 

-0 6016 

-1 2032 
(1 0946x10-^)* 

Cir 

0 1893 

0 1867 

0 1893 

(2.6677x10’^) 

Cip 

-0 1444 

-0 1108 

-0 1444 
(1 1420x10-^) 

Cl8r 

0 0318 

0 0363 

0 0318 
(1 3400x10'^) 

C„p 

-0 2101 

-0.2347 

-0 2101 
(1 4870x10^) 

Cnr 

-0 8050 

-1.1188 

-0 8050 
(3 5481x10''^) 

Cnp 

0 1927 

0 2189 

0.1927 

(1.5513xl0‘'^) 

Cn6r 

-0 1649 

-0.1405 

-0.1649 

(1.8218x10-*) 

Cyp 

0 9903 

1.0771 


Cyr 

1 5937 

1.2233 


CyP 

-0.7726 

-0.7691 


Cy5r 

0 3016 

0 2983 

0.3016 

(8.6147x10-'') 


*Cramer Rao Bounds 









Table 3.9 Parameter Estimates from Raw Flight Data using Delta Method of FFNN 
(Lateral-Directional Motion) 


Parameter 

Estimates 


FLT5/TP1 5 
(V, = 160 KIAS) 

FLT5/TP9 14 
(V.= 180 KIAS) 

FLT5/TP17 22 
(V, = 200 KIAS) 

FLT5/TP24 28 
(Vi = 220 KIAS) 

Cip 

-1.1646 
(0 1832)* 

-1.1585 
(0 2845) 

-1 1935 
(0 0856) 

-0.9354 

(0.1898) 

Cir 

0 4269 
(0 2550) 

0 1143 
(0 3794) 

0 3841 
(0 0507) 

0.3484 

(0.0834) 

Cip 

-0 1658 
(0 0445) 

-0.1132 
(0 0800) 

-0 1873 
(0 0346) 

-0.1285 

(0.0206) 

Cl6a 

-0 4542 
(0 0879) 

-0.3236 

(0.0875) 

-0.3473 
(0 0460) 

-0 3302 
(0.0437) 

Cl6r 

0.0227 
(0 0430) 

0.0564 
(0 0910) 

0 0339 
(0.0088) 

0.0307 
(0 0083) 

Cnp 

-0 2481 
(0 0908) 

-0.1812 
(0 0971) 



C„r 


-0 7074 
(0.0460) 

-0.6328 

(0.0828) 

-0.5505 

(0.2022) 

C„p 

0.2376 
(0 0408) 

0.2588 

(0.0263) 

0 2173 
(0.0263) 


Cn6a 

-0.0730 

(0.0293) 

-0.0509 

(0.0361) 

-0.0343 

(0.0244) 


Cn5r 

-0 1293 
(0.0423) 


-0.1330 

(0.0077) 

QQn|S9|||H 

Cyp 

3.0455 

(1.3590) 

0.7622 

(0.2334) 

1.0884 

(0.3803) 

0.3726 

(0.1989) 

Cyr 



BmHH 

2.0186 

(0.5459) 

CyP 

-0.8169 

(0.0919) 

-1.0063 

(0.0756) 

-0.9227 

(0.0835) 

-0.9406 

(0.0580) 

Cy5a 

0.8754 

(0.3232) 

0.3081 

(0.0890) 

0.2527 

(0.1015) 


Cy5r 

0.3633 

(0.1249) 

0.1314 

(0.1473) 

0.3144 
(0 0789) 

0.3161 

(0.0892) 


* Standard Deviation 
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Table 3.10 Estimates of Parameters Cio , Cno and Cyo for Raw Flight Data 
(Lateral and Directional Motion) 


Parameter 

Test Points 

FLT5/TP1 5 
(V,= 160KIAS) 

FLT5/TP9 14 
(Vi = 180 KIAS) 

FLT5/TP17 22 
(Vi = 200 KIAS ) 

FLT5/TP24 28 
(Vi = 220 KIAS) 

C,o 

-0.0086 

-0.0100 

-0.0136 

-0.0079 

Cno 

0.0114 

0.0154 

0.0130 

0.0132 

Cvo 

-0.0631 

-0.0788 

-0.0723 

-0.0722 
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Table 3.11 Parameter Estimates from Corrected Flight Data using Delta Method of 
FFNN (Lateral-Directional Motion) 


Parameter 

Estimates 


FLT5/TP1 5 
(Vi = 160 KIAS) 

FLT5/TP9 14 
(Vi = 180 KIAS) 

FLT5/TP17 22 
(V, = 200 KIAS) 

FLT5/TP24 28 
(Vi = 220 KIAS) 

Cip 

-1 2032 
(0 1925)* 

-0.9370 

(0.0858) 

-1 0820 
(0 0802) 

-1.0928 

(0.0746) 

Clr 

0 1893 
(0 0717) 

0.0648 
(0 2318) 

0 2703 
(0 1889) 

0 2497 
(0.1578) 

Cip 

-0 1444 
(0.0559) 

-0 1214 
(0.0658) 

-0 1299 
(0 0271) 

-0.1137 

(0.0109) 

Cl8a 




-0.3348 

(0.0634) 

Cl5r 

0 0318 
(0.0295) 

0.0497 

(0.0626) 

0 0367 
(0 0066) 

RQ^|||||||H 

Cnp 

-0 2101 
(0 0586) 

-0 2618 
(0.0414) 

-0.2290 
(0 0654) 


Cnr 

-0 8050 
(0.0975) 

-0.7993 

(0.0898) 

-0.6380 
(0 1254) 

loQiBm 

C„p 

0 1927 
(0 0341) 

0.2417 
(0 0374) 

01858 
(0 0454) 


Cn6a 

-0.0374 

(0.0457) 

-0.0621 

(0.0200) 

-0.0614 

(0.0304) 

Q^n|||l|[[H 

IE9BSnHH 

Cn5r 

-0.1649 

(0.0206) 

-0.1608 

(0.0537) 



Cyp 

|^^^^■|||||■ 



|QQESgHH||||||| 

Cyr 





Cyp 

-0.7726 

(0.0866) 

-0.8529 

(0.0804) 

-0 7864 
(0 0961) 

-0.8064 

(0.0809) 

Cy6a 

0.5976 

(0.1816) 

0.4375 

(0.1415) 


0.2461 

(0.1095) 

Cy5f 

0.3016 

(0.2391) 

0.2981 

(0.1655) 

0 3616 
(0.1097) 

0.4173 

(0.1319) 


* Standard Deviation 
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Table 3.12 Estimates of Parameters Q© , Cn© and Cy© for The Corrected Flight 
Data (Lateral-Directional Motion) 


Parameter 

Test Points 

FLT5/TP1 5 
(V,= 160 KIAS) 

FLT5/TP9 14 
(V.= 180 KIAS) 

FLT5/rP17 22 
(V, = 200 BQAS ) 

FLT5/TP24 28 
(Vi = 220 KIAS) 

C,o 

-0.0020 

-0.0039 

-6.3368x1 0"^ 

0.0014 

Cno 

9.2389x10'^ 

7.3766x10"^ 

-0.0019 

-0.0017 

Cyo 

0.1233 

0.0935 

0.0930 

0.0835 




CHAPTER 4 


CONCLUSIONS AND RECOMMENDATIONS 

4.1 CONCLUSIONS 

In the present work, the Delta method using FFNNs has been used to estimate the 
stability and control derivatives from the real flight data, pertaining to longitudinal and 
lateral-directional dynamics. The Delta method has been used to extract the parameters 
from the raw as well as the corrected flight data. The results of the parameter estimates 
obtained from the Delta method in each case, show good consistency with low standard 
deviation. Further, a close match is obtained between the estimated and actual 
aerodynamic coeflhcients as well as between the estimated and measured responses. This 
implies that the Delta method estimates are fairly accurate. The Delta method has been 
found to be a simple and straight forward method that can be set up very quickly. 
Moreover, unlike the conventional output error methods (ML), the Delta method does not 
require any knowledge of the initial guess values of the parameters. Further, the output 
error methods will require integration of system equations of motion. This step is 
bypassed in the Delta method. This is a distinct advantage of the Delta method. In foct, 
it can be used to give a quick look at the values of the parameters, if so be desired. 
Subsequently, with little additional fine tuning of the network, reasonably good estimates 
can be obtained, as seen in the present work. This method has been proven to be a robust 
method, which can extract parameters even from the raw flight data usmg the measured 
values of the motion variables that are not corrected for bias and scale factor errors. 
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4.2 RECOMMENDATIONS 


Based on the experience gathered during the present work, a few 
recommendations are suggested for foture work: 

1. For the parameter estimation for lateral-directional motion, instead of giving 
separate 5a and 5r inputs, a combined input should be tested, wherein the 5r input is 
given first and following 1-2 second delay the 5a input should be given. 

2. The data analysis in the present work for the Delta method (on real flight data) 
was done only in the time domain. It would be of great interest, to see how the results 
compare, with the same analysis being done in the frequency domain. 

3. The different architectures of FFNN in the present work were obtained by using a 
few functions/features of the Matlab Neural Network Toolbox. This aspect can be 
investigated in greater detail by using all the different functions to see the improvements 
if any, in the results. 

4. The criterion used for freezing the optimal FFNN architecture during training 
phase was the minimum tr ainin g MSB criterion. Other criteria could be investigated 
which may possibly lead to better training and hence improved prediction c^>abi]fry of 
the FFNN. 

5. Future projects, along the lines of the present work, could be jointfy undrataken 
by IIT Kanpur and ASTE, AF, Bangalore to explore the high angle of attack regime, as 
well as for parameter estimation of open loop xmstable airaaft. 


Ill 
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APPENDIX ’A* 


FEED FORWARD NEURAL NETWORK 

Artificial Neural Network can be described as an attempt by researchers to mimic 
the functioning of the human brain Meinsky and Papert^^ carried out a rigorous analysis 
of the perceptron and pointed out certain limitations of neuron models, the publication of 
these results brought the research in this area to a halt. Later, the basic idea of back 
propagation propounded by Werbos^^, and work of Hopfield^^ revived interest in neural 
networks. Presently, research on neural network is being actively pursued in many 
disciplines, ranging fi’om neurobiology and physiology to engineering sciences 

An artificial neural network (ANN) has a parallel and distributed processing 
structure and consists of processing elements called neurons. These neurons are 
interconnected with weighted unidirectional connections. Each neuron is composed of 
multiple inputs, one output, local memory and activation function. The inputs carry the 
weighted output of the directly connected neurons. The incoming information of a 
neuron is processed by the associated activation function (such as a sigmoidal function). 
The output is then distributed to other neurons as inputs. A variety of ANN 
configurations has been developed. One of the most popular ANN configurations is the 
feed forward neural network (FFNN) shown in Fig A.I. 

Typically, the FFNN is made up of source nodes that constitute the input layer, 
and one or more hidden layer, and an output layer. The input signal gets propagated 
through the network in a forward direction, on a layer by layer basis. Each node in the 
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Fig. A. I 


Schematic Representation of Feed Forward Neural Network 
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layer (input or hidden) is connected to each node in the neM layer (either the hidden or 
the output) through the connective weight Every node in the input layer and the hidden 
layer is biased 

A Feed Forward Neural Network has three distinct characteristics. 

i) Each neuron perfornis aggregation of its input and this aggregated information is 
passed through the nonlinear activation function One of the most commonly used 
activation function, the sigmoidal (logistic) function is defined as 

V. 1 

* 1 + exp(-Xj IX) ^ 

where x, is the net internal activity level of neuron i, Y, is the output of the neuron and X 
is the logistic gain function 

ii) The network contains one or more layers of hidden neurons that are not part of 
either the input or the output layer. These are not part of either the input or the output 
layer. These hidden neurons are primarily responsible for learning complex phenomena 
by progressively extracting more meaningful characteristics fi’om the input pattern. 

iii) The neurons in the architecture are progressively connected to the neurons of the 
next layers through connective weights. Although the model of the system to be 
approximated is, to a large extent, encoded in these weights, there is no direct 
relationship between the weights and parameters to the system to be modeled. 

FFNNs have been applied successfully to model complex physical phenomena by 
training them in a supervised manner with a popular algorithm known as the back- 
propagation algorithm (BPA)^^. The formulation of back-propagation algorithm provided 
a computationally efficient method for training of FFNN. Fundamentally, the BPA 
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consists of a forward and a backward pass through the different layers of the network. 
During the forward pass, the input vector is applied to the nodes of the network and its 
effect propagates through the network layer by layer The connective weights are all kept 
fixed during the forward pass On the other hand, in the backward pass, the weights are 
updated in accordance with the error correction rule Specifically, the actual response of 
the network is subtracted from the desired response to produce an error signal This error 
signal is then propagated backward against the direction of the connective weight-hence 
the name 'back-propagation algorithm'. 

Back-propagation is a specific technique for implementing gradient descent in 
weight space for feed forward neural network In a practical application of back- 
propagation algorithm, learning is achieved from the many presentations of a specific 
training example to the FFNN. For a given training set, back-propagation learning may 
proceed in Pattern mode or Batch mode In pattern mode of back-propagation learning, 
the network weights are updated sequentially as the training data set is being presented 
Then the net weight change, averaged over the training set of pattern, is used to update 
the network weights. On the other hand, in the batch mode, weight updating is performed 
after presentation of all the training examples. 

It may be pointed out that any back-propagation training algorithm is recursive in 
nature and it needs repetitive training sessions to achieve the required learning. There are 
many network parameters (also called the influence parameters or the tuning parameters) 
like the learning rate, the momentum constant, the number of hidden layers and the 
number of neurons in each of the hidden layers, number of iterations for which the 
training data is passed through the network, the logistic gain constant of the sigmoidal 
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function, etc that effect the accuracy of functional mapping between the input and the 
output variables 

The most general network performance used is the "Mean Square Error (MSE)" 
function MSE measure the network's performance according to the mean of squared 
errors Thus MSE can be defined as 

MSE = E Z [yi(j) - (A 2) 

m * n j=i i=i 

where y and x are respectively the desired (known) and the computed (predicted) outputs 
of the neural network, n is the number of data points, and m is the number of the output 
variables. If only one output variable is to be trained, then MSE is defined with m=l 
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APPENDIX B 


MAXIMUM LIKELIHOOD METHOD (MLE) 

In order to provide a comparative study for the efficacy and reliability of the 
FFNN method, the longitudinal parameters were also estimated using the conventional 
'Maximum Likelihood Method' The programs were written in Matlab using the MLE 
algorithm The algorithm was derived as follows 

Let a mathematical model of a system be represented by the following equations 
X(ti) = AX(t,) + BU(t.) (B.l) 

Z(ti ) = CX (ti) + Gm (B 2) 

where X is the state vector, U the control vector, Z the measured (observation) response 
and nj is the noise vector respectively. The matrices A and B contain the unknown 
stability and control derivatives (parameters). If there is no state noise and the matrix G 
is known, then the maximum likelihood estimator minimizes the cost function, 

m - i 2[Z(‘i) - Ze(t,)f (GG'^r‘[Z(t,) - ZsOi)] (B.3) 

2 1=1 

where GG^ is the measurement noise covariance matrix and Z 0 (tj) is the computed 

response estimate of Z at h for a given value of the unknown parameter vector 0. The 
cost function is a function of the difference between the measured and computed time 
histories In our case, we have taken GG^ = I for the cost function. To minimize the cost 


122 



function J(0), we can apply the Newton-Raphson algorithm, which chooses successive 
estimates of the vector G of unknown coefficients If k is the iteration number, then 

0k+l = 9k -[''i^0j(Qk)r^[VeJ(0k)] (B 4) 

where the first gradient is defined by 

VeJ(e) = - z[Z(ti) -Ze(ti)f (GG'r)-i[VeZe(ti)] (B 5) 

1=1 

and the second gradient is approximated by 

V^eJO) = I [VeZeCtof (GGTr‘[VeZe(t,)] (B 6) 

1=1 

The matrix GG^ is approximated by a diagonal matrix, with the diagonal elements given 
by, 

° [Z(t,)-ZB(ti)]Z(ti)-Ze(ti)f 
i=i N 

where N is total number of data points. The term VeZeCtj) occurring in the first and the 

second gradients is evaluated by finite difference approximation for each of the 
sensitivity coefficients present in it. For example, to evaluate sensitivity coefficient 

, the equations of motion are solved to obtain a for Cmq + ACmq and Cmq - ACmq . 

^mq 

Let the respective values of a be and a , then we have 


da 

^mq 


a^ -a 

2ACjiiq 
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4 i Ilf 


The reliability of the estimated parameter in terms of the Cramer-Rao bounds is given by 

the square roots of the diagonal elements of [V0j(0]s.)]''\ which is any way evaluated as 

part of the ML algorithm In the present work, for estimating the longitudinal 
parameters, the vector of unknown parameters is given by, 

Q ~ [CLa>^Lq!^L5e5^ma’^mq3^m6e] 

with the Z(ti) and Z0(tj) responses being a and qc/ 2 u, while the control vector has 5 e as 
the input 
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